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Abstract 


White dwarfs are a class of stars with unique physical properties. They 
present many challenging problems whose solution requires the application of 
advanced theories of dense matter, state-of-the-art experimental techniques, 
and extensive computing efforts. New ground- and space-based observatories 
will soon provide an increasingly detailed view of white dwarf stars and reveal 
new phenomena that will challenge our models. This review is an introduc- 
tion for researchers who are not in the field of white dwarf astrophysics with 
the intent to entice them to contribute their expertise to advance our knowl- 
edge of these exotic stars. We discuss a wide variety of currently unsolved 
or partially resolved problems that are broadly related to equations of state, 
transport processes and opacities. 
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1. Introduction 


At the end of their long lives, the vast majority of stars become white 
dwarf stars, and such is the fate of our Sun. These common denizens of 
our Galaxy have been studied by astrophysicists for about a century. White 
dwarfs are the sites of exotic physical conditions that are not usually encoun- 
tered in normal stars. From the tenuous surface with a temperature of the 
order of 1 eV to the center where it reaches ~ 100 eV and the density rises to 
~ 10° ¢cm~?, white dwarfs present a number of challenging physics problems 
that need to be addressed to properly model their structure, evolution and 
spectra. 
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Historically, theoretical work in dense plasma physics has found fertile 
applications in white dwarf models. A famous example is one of the first ap- 
plications of quantum mechanics to astrophysics—by R. H. Fowler in 1926— 
which revealed that the abnormally small radius of white dwarfs follows from 
the equation of state of the degenerate electron gas. Further progress in un- 
derstanding white dwarfs required the development of models for equations 
of state, phase diagrams, radiative opacities, electron conductivity and diffu- 
sion coefficients in complex physical regimes. Some of these models remain 
approximate and sometimes inadequate to account for the observations. On 
the other hand, the composition of white dwarf matter at a given depth in the 
star is typically mono-elemental or a binary mixture, which offers a degree 
of simplification compared to other types of stars. 

Thousands of white dwarfs have been studied in detail and their ranks 
have grown considerably to about a quarter million, most of them identified 
in the last few years. The astrophysics of white dwarfs is a mature field 
and there is an abundance of exquisite data that calls for new and refined 
physics models and also provides actual constraints on theory. The spectra 
of white dwarfs reveal conditions at the surface, where the gas can be dense 
enough to be interacting, with observable effects on the chemical/ionization 
equilibrium, the equation of state, and the continuum and line opacities. 
The pattern of surface abundances reveals the crucial importance of inter- 
species diffusion in the star’s high gravity field. Some white dwarfs, like 
many other types of stars, display periodic variations in their brightness that 
arise from global oscillations. Most white dwarfs go through at least one 
pulsation phase during their evolution. In a manner entirely analogous to 
seismology for the Earth, the study of these global oscillation modes—known 
as asteroseismology—probes the interior structure of the white dwarfs in 
remarkable detail. 

Beyond their intrinsic astrophysical interest, white dwarfs ought to com- 
mand the attention of physicists. They present exotic physical conditions 
that are not seen on Earth and that challenge theories of dense matter and 
even not-so-dense matter. While state-of-the-art experimental techniques 
are nibbling at the edges of the least understood physical regimes, there is 
an abundance of astrophysical data that bears strongly on the underlying 
physical theories and models. The current quantity and precision of the 
astrophysical data are expected to rise dramatically in the next few years, 
thanks to new space-based and ground-based observatories and large-scale 
surveys of the sky. The field is ripe with interesting problems whose solution 


can be validated against data. 

This review is aimed at scientists with an interest in modeling and mea- 
suring the properties of materials under unusual conditions. Enough infor- 
mation is given about white dwarfs to provide context and motivation, but 
the emphasis is on physical processes and unsolved problems of immediate 
interest. We introduce basic concepts of stellar structure, evolution and at- 
mospheres as needed for clarity and context but these are kept to a minimum. 
Several excellent monographs at the graduate study level provide introduc- 
tions to stellar astrophysics and we cite recent reviews of white dwarf stars 
that expound on this rich field of research. 

The review is structured in layers, with each main section delving deeper 
into the phenomenology and physics of white dwarfs. We first describe the 
astrophysical context and significance of white dwarfs in broad terms (Sec- 
tion 2). This is a high-level introduction to white dwarfs as a class of stars. 
The principal physical processes responsible for their structure, evolution 
and spectral appearance are introduced in Section 3. This is followed by a 
discussion of the physical regimes encountered in white dwarfs that sets the 
context for the constitutive physics needed to construct white dwarf mod- 
els (Section 4). The stage is then set for introducing several contemporary 
problems in the physics of white dwarfs (Section 5), some whose current treat- 
ment needs refinement and some that remain unsolved. This section covers 
the more detailed physical theory of many topics including phase diagrams 
of dense plasmas, transport coefficients in partially ionized dense plasmas, 
collision-induced absorption opacity, line profile theory, and simulations of 
convection. Finally, we conclude with a list of problems that we hope will 
compel readers to contribute to the advancement of the astrophysics of white 
dwarf stars (Section 6). 

Many phenomena involving white dwarfs, particularly in interacting bi- 
nary star systems, are quite complex and are part of exciting branches of 
modern astrophysics. This review focuses on the somewhat simpler case of 
single white dwarfs as the potential for gains in understanding by the appli- 
cation of physical theory and models is, in our opinion, greatest and more 
immediate. 


2. White Dwarfs as Stars 


White dwarfs represent the final evolutionary stage of stars less massive 
than ~ 10M, ' —more massive objects will collapse into neutron stars or 
black holes. This includes about 97% of stars in the Milky Way. For most 
stars, this ultimate state is reached after fusing H into He (first on the main 
sequence and later on the red giant branch), then fusing He into C and O (first 
on the horizontal branch and later on the asymptotic giant branch, AGB), 
and finally losing their outer layers until only a dense C-O core remains. 
Without any nuclear energy source, except for residual H burning in some 
cases (Section 3.1), these dense stellar embers simply cool down for the rest 
of time. 


2.1. Anatomy of a white dwarf 


A typical white dwarf has a mass of 0.6 Mo and a 0.013 Ro (= 1.4 Re) 
radius, implying an extremely high density (9 ~ 10° gcm~? throughout most 
of its structure). Under those conditions, matter is fully ionized and the 
electrons are completely degenerate throughout the bulk of the star. The 
pressure of this degenerate electron gas dominates the structure of the star 
and prevents it from collapsing under the influence of its strong gravity’. The 
hydrostatic equilibrium between the inward pull of gravity and the outward 
push of the electron degeneracy pressure leads to the well-known mass—radius 
relation of white dwarfs. For a given chemical composition, the radius of a 
white dwarf is uniquely determined by its mass and temperature. Since the 
degeneracy pressure increases rapidly with density, the higher the mass of 
a white dwarf, the smaller its radius. For white dwarfs supported by the 
pressure of fully degenerate but non-relativistic electrons, the mass-radius 
relation is M, o R;z?, where M, and R, are the mass and radius of the 
star, respectively. There is however a limit to the mass that electron degen- 
eracy pressure can support. At very high densities, the electron gas becomes 
relativistic and the equation of state softens (i.e., goes from P « ne? in 


the non-relativistic case to P « n2/* in the strongly relativistic limit, where 


‘Convenient astronomical units are the mass, radius and luminosity of the Sun, Mo = 
1.988 x 10°° kg, Ro = 6.957 x 10° m and Lo = 3.828 x 107° W, respectively, and the radius 
of the Earth, Re = 6.378 x 10° m. 

In contrast, the lower-density main sequence stars are supported by thermal pressure, 
a negligible quantity in white dwarf interiors. 


P is the pressure and n, the electron density). This reduces the capacity of 
electron degeneracy pressure to support the white dwarf, leading to an upper 
limit on the mass of a stable white dwarf. This limit of = 1.4 Mo is known 
as the Chandrasekhar mass [1]. 

For a typical white dwarf, the C and O leftover after the post-main se- 
quence phases of the evolution form a core that represents 99% of the total 
mass (Figure 1). Owing to the large thermal conductivity of its degenerate 
electrons, this core is nearly isothermal, with T = 10? to 10*+eV depending 
on the age of the white dwarf. The remaining outer 1% is made of residual 
He and H that was not fused during the earlier evolutionary phases. An He 
envelope (~1% of the total mass) surrounds the C-O core, which is itself sur- 
rounded by a thin H layer (< 0.01% of the total mass). The exact thickness 
of those H and He layers depends on details of nuclear burning during the 
late phases of stellar evolution. The uppermost layers of the H/He envelope 
form the atmosphere of the white dwarf, the only region we can directly ob- 
serve. Its temperature ranges from T = 10 to 0.1 eV, and its density remains 
<1gcem7?. While the H/He/C-O structure outlined above is representative 
of most white dwarfs, there are many exceptions. Notably, about 20-25% of 
white dwarfs have little or no H present [2] because they burned virtually all 
their H during previous evolutionary phases [3]. 

The highly stratified chemical structure of white dwarfs is due to their 
compact nature. Their high densities confer on them a very strong surface 
gravity (g ~ 10° cms~?, compared to 3 x 10*cms~? for the Sun), with the 
consequence that heavier elements quickly sink toward the core and lighter 
ones float to the surface [4]. While gravitational settling is the dominant force 
shaping the chemical structure of a white dwarf, it is often in competition 
with other mechanisms such as convection [5, 6, 7, 8], accretion [9, 10], and 
radiative levitation [11]. In any case, the distinct H/He/C-O layering of 
white dwarfs greatly simplifies their modeling compared to most other types 
of stars, where difficult-to-model convective mixing usually plays a large role. 

A C-O core composition is not a universal property of white dwarfs. If 
the progenitor star is close to the 10 Mo limit, it may reach high enough 
temperatures to fuse C, thereby leading to an O-Ne white dwarf. The so- 
called ultra-massive white dwarfs (M, > 1.1 Mo) are thought to have O-Ne 
cores [12] (although recent evidence points to the existence of ultra-massive 
C-O white dwarfs as well [13, 14, 15, 16, 17]). At the other extreme, very 
low-mass progenitors (< 0.6M.) will not reach the conditions required to 
fuse He into C and O, leading to a He-core white dwarf. In isolation, the 
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Figure 1: Chemical layering in a representative model of a 0.6 Me H-dominated at- 
mosphere (DA) white dwarf with a surface temperature Tig ~ 12000K. The X-axis 
uses a mass coordinate where m(r) is the mass within a sphere of radius r, hence 
log q = log(1 — m(r)/M,) is the log of the mass fraction outside of radius r. The center of 
the star is at the left, the surface is towards the right. The atmospheric layers, located at 
log g < —15, are not shown here. 
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Figure 2: Mass distribution of white dwarfs, based on the sample of Ref. [21]. 


low-mass progenitors of He-core white dwarfs are expected to remain in the 
H-fusing main sequence phase longer than the age of the Universe. Any He- 
core white dwarf that we observe today likely formed by mass loss in a close 
binary system [18, 19]. This is notably the case of the extremely low-mass 
(ELM) white dwarfs (MM, < 0.3. Mo). While white dwarfs with M, + 0.2 Mo 
all the way to 1.35 Ms have been identified, the vast majority have masses 
remarkably close to 0.6 Mo. Indeed, ~ 60% of white dwarfs have a mass 
comprised between 0.5 and 0.7 Me [20, 21] (see Figure 2), and + 95% are 
between 0.45 and 1.1 Mo, where a C-O core is expected. 

Around 10% of known single white dwarfs are found to harbor a detectable 
magnetic field. The strength of this field spans many orders of magnitudes, 
from 1kG (= 0.1T) for the weakest fields detected to 1000 MG for the most 
strongly magnetic white dwarfs known (for comparison, sunspots can reach 
3 kG) [22]. The origin of those fields remains a fundamental unsolved problem 
of stellar astrophysics. A remnant fossil field inherited from their progenitors, 
the generation of a field during close binary evolution and stellar mergers, 
and a dynamo process in the white dwarf itself are all possible explanations. 


2.2. White dwarf cooling 


Once a star reaches the white dwarf state, it simply cools down for the rest 
of its evolution. Accurately modeling this cooling process has been one of the 
central goals of white dwarf research for over half a century [23, 24, 25, 26, 27]. 
While the evolution of white dwarfs is relatively simple compared to other 
astrophysical objects, it demands a precise description of a variety physical 
processes, which take place under extreme conditions seldom found elsewhere 
in the Universe. The rate at which a white dwarf cools down depends on the 
neutrino production rate in its core, the thermal conductivity of its core and 
envelope, the radiative opacity of its atmosphere, element transport in its 
interior, and the physics of phase transitions. The problem of white dwarf 
cooling is therefore far from trivial, but the unique potential of white dwarf 
cosmochronology (Section 2.5) motivates theoretical and experimental efforts 
to improve our description of the physical processes at play. 


2.8. White dwarf atmospheres 


All the information we receive from white dwarfs comes from their at- 
mospheres, making this very thin surface layer disproportionately important 
to our overall understanding. The atmosphere also plays the role of outer 
boundary condition for interior and evolution models, providing further mo- 
tivation to achieve an accurate description of this region [25, 28, 29]. 

A first way to characterize the atmosphere of a white dwarf is through 
photometry which consists of measuring the spectral flux density f, emitted 
by the star and received on Earth within a given bandpass. This emergent 
flux mostly comes from a region known as the photosphere, where 50% of the 
photons escapes from the atmosphere without being absorbed or scattered. 
By measuring the flux within a series of different bandpasses, the spectral 
energy distribution of the star is constrained. This spectral energy distribu- 
tion is used to determine the surface temperature (reported as the effective 
temperature Ty, which corresponds to the temperature of a blackbody that 
would emit the same total amount of radiation). In first approximation, the 
spectral energy distribution is given by the spectrum of a blackbody with 
temperature Tx¢. However, strong deviations from this approximation are 
common due to the various opacity sources in the atmosphere, and in prac- 
tice detailed atmosphere models are used [30]. In some cases, the atmospheric 
composition (normally dominated by H or He) can be constrained from the 
spectral energy distribution alone [31, 32, 33, 34]. 


If the distance D of the white dwarf to the Earth is known (which is now 
the case for most known white dwarfs, thanks to the Gaia space observatory), 
the mass of the white dwarf can also be inferred. The radius R, of the white 
dwarf is obtained using 


2 
pobscrved _ TRY intrinsic ( 1 ) 
Xr ~~ D 2 aN ? 


where fi sic is the intrinsic flux given by a model atmosphere, and fesrved 


is the flux measured at the Earth. From R,, the mass of a white dwarf (and 
surface gravity, g = GM,/R?) is then directly obtained from the theoretical 
mass-—radius relation. 

Another way through which the atmosphere of a white dwarf can be char- 
acterized is with spectroscopy. Spectroscopy generally enables a more precise 
characterization of the atmosphere as it unlocks the wealth of information 
contained in spectral lines and molecular bands. It gives access to the de- 
tailed chemical makeup of the atmosphere, to its temperature and pressure 
conditions (through the strengths and widths of the different spectral fea- 
tures), and to its magnetic field (via the Zeeman effect). The analysis of 
those spectral lines and molecular bands requires a detailed description of 
a variety of phenomena that affect the spectral line shapes (Stark effect, 
pressure broadening, etc.), as well as a good model for the vertical stratifica- 
tion of the atmosphere since the spectrum is formed from the contributions 
of many distinct layers that each have their own temperature and pressure 
conditions. 

The diversity of atmospheric compositions and surface temperatures gives 
rise to a variety of “spectral types”. The most important ones are listed in 
Table 1, and an example of each one is shown in Figure 3. The DA spectral 
type is by far the most common and is characterized by H lines from the 
Balmer series in the visible portion of the electromagnetic spectrum. Then, 
there are the types DB and DO, which show He and He? spectral lines, 
respectively. DB and DO white dwarfs have very little or no residual H, 
which explains the visibility of He/He* lines. At low temperatures (below 
Teg * 11000K for He and 5000K for H), He and H spectral lines are no 
longer detectable in the visible portion of the spectrum, yielding a continuous, 
featureless spectrum. Those objects are known as type DC. There are also 
objects with C features in their spectra, known as DQ white dwarfs. Those 
C features can either be C/Ct spectral lines (above Tog = 10000 K) or C2 


Spectral type Defining characteristics Origin 


DA H spectral lines Standard single-star evolution 

DB He spectral lines H deficient 

DO Het spectral lines H deficient 

DC No spectroscopic features Cool atmosphere 

DQ C or Ct spectral lines or C2 ~—-C dredge-up or stellar merger 
molecular bands 

DZ Spectral lines from elements Accretion of rocky debris 
heavier than C (e.g., Ca, 
Mg, Fe) 


Table 1: The most common white dwarf spectral types. 


molecular bands (at lower temperatures). How can C be found in white dwarf 
atmospheres if gravitational settling is so efficient? Two distinct evolutionary 
channels are thought to explain the DQ phenomenon. On the one hand, hot 
and massive white dwarfs with C-dominated atmospheres (known as hot DQs 
[35]) are thought to be formed after the merger of two white dwarfs in a binary 
system [36, 37]. On the other hand, the presence of C in cooler objects with 
C-polluted (but He-dominated) atmospheres is explained by the dredge-up of 
C from the deep interior by a convection zone in the He layer [6, 38]. Finally, 
there are the DZ white dwarfs, characterized by spectral lines from elements 
heavier than C. In an exciting recent development, it is now accepted that 
the presence of low-Z elements (e.g., Ca, Mg, Fe) in their atmospheres is 
due to the recent accretion of rocky debris from their own planetary system 
[10, 39, 40]. This scenario is supported by the identification of debris disks 
around metal-polluted white dwarfs [41, 42, 43], many of which have excess 
infrared emission due to warm dust [39], the detection of recurring transit 
events caused by planetary debris [44, 45, 46, 47] and the measurement of 
X-ray emission induced by the accretion event [48]. 


The spectral type of a white dwarf is not a direct indication of its atmo- 
spheric composition. For example, a DA white dwarf can have a substantial 
amount of He in its atmosphere and still show only H lines, and most DQ and 
DZ white dwarfs have atmospheres dominated by He. The spectral type only 
denotes which species dominates the spectrum. Note also that the spectral 
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Figure 3: Example spectra of DA, DB, DO, DC, DQ, and DZ white dwarfs. All spectra 
are normalized at 5000 A and vertically shifted for clarity. Spectra taken from [49, 50, 51, 


52, 2]. 
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type is not a fixed property of a white dwarf. Not only does the contin- 
uously decreasing temperature of the atmosphere determine which spectral 
features can be detected, but, more interestingly, the chemical composition 
of the atmosphere also changes. A He-dominated atmosphere can transform 
into a H-dominated atmosphere, and vice versa, due to transport mecha- 
nisms (chiefly, convection) that compete with the ever-present gravitational 
settling. This constitutes a unique window into processes that occur below 
the directly observable portion of the atmosphere, and it is the subject of 
intense observational and theoretical efforts [8, 32, 37, 53, 54, 55, 2, 33, 34]. 


2.4. Pulsating white dwarfs 


Little can be learned about the deep interiors of white dwarfs from the 
analysis of their atmospheres alone. Fortunately, the pulsations of white 
dwarfs can be used to infer the properties of their deeper layers. In white 
dwarfs, the pulsations are gravity modes where buoyancy is the restoring 
force, with periods of © 100—2500s. Each pulsating star vibrates in a unique 
set of normal modes, the frequencies of which are determined by its chemi- 
cal stratification, rotation profile, mass, radius, temperature, and magnetic 
field. Asteroseismology is the study of those normal modes with the aim 
of extracting information on the inner structure of stars. Excellent reviews 
are devoted to white dwarf asteroseismology [56, 57, 27, 58]; here we only 
provide the context needed for the following sections. 

Only about ~ 2 — 3% of white dwarfs pulsate, but most will pulsate at 
some point in their evolution. The cause of the onset of pulsations in white 
dwarfs is the partial ionization of the main constituents of their envelopes. 
Partial ionization leads to a formidable increase of the envelope opacity, which 
impedes the outgoing radiation flow and causes pulsational instabilities. The 
composition of the envelope determines when the partial ionization regime 
is reached. Hence, DB white dwarfs pulsate between Tog = 32000K and 
22000 K due to the partial ionization of He (forming the DBV pulsator class), 
and DA white dwarfs pulsate in the Teg = 12500 — 10500 K range when H 
is partially ionized (forming the DAV class). There are at least six different 
types of white dwarfs pulsators currently known, each falling within a narrow 
range of parameters—predominantly T.4¢—that defines an instability strip. 
All of those pulsator classes correspond to different phases of the evolution, 
thereby offering unique peeks at the inner structure of white dwarfs through 
specific windows during their evolution. 
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2.5. White dwarfs and modern astrophysics 


Besides being interesting objects in themselves, white dwarfs are also use- 
ful to other branches of astrophysics. As they age, the observable properties 
of the H-burning main sequence stars remain almost unchanged, rendering 
age determinations very uncertain [59]. In contrast, because of their contin- 
uous cooling, the observable properties of white dwarfs change drastically as 
they age, making them uniquely precise cosmic clocks. They have been used 
to measure the age of a variety of stellar populations [60, 61, 62, 63, 64, 65, 66] 
and to reconstruct the history of our Galaxy |67, 68, 69, 70]. To determine 
the age of a white dwarf, one must first determine its mass, temperature, and 
atmospheric composition, which requires atmosphere models. Those atmo- 
spheric parameters are then compared to the predictions of cooling models to 
obtain a white dwarf cooling age (to which the appropriate pre-white dwarf 
lifetime can be added to obtain the total stellar age). Precision white dwarf 
cosmochronology therefore requires accurate atmosphere and cooling models. 
Considering the current limitations of those models, the uncertainty on the 
ages of very old white dwarfs (~ 10 Gyr?) likely exceeds 1 Gyr. Improving 
this figure is one of the chief goals of white dwarf astrophysics. Achieving 
an accuracy of © 4% would provide useful constraints to distinguish between 
competing cosmological models [71]. 

As a fossil record of previous stellar evolutionary phases, the interiors of 
white dwarfs contain precious information on the physical processes at play 
in other types of stars. As such, white dwarf asteroseismology is a promising 
tool to answer current questions in other subfields of stellar astrophysics. 
In particular, the C-O abundance profiles of white dwarfs contain critical 
information on the difficult-to-model convective mixing at the boundary of 
the He-fusing cores of red giant stars [72, 73, 74], as well as on the poorly 
constrained !?C(a,~)!®O reaction rate, dubbed “the holy grail of nuclear 
astrophysics” |75, 76, 77, 78, 79, 80]. Furthermore, the internal structure 
of high-mass white dwarfs that descend from double white dwarf mergers 
(including hot DQs) [81] may provide valuable constraints for numerical sim- 
ulations of such mergers. Double white dwarf mergers are one of the channels 
through which a type Ia supernova (the powerful explosion of a C-O white 
dwarf that exceeds the Chandrasekhar limit) can occur [82, 83]. High-mass 
white dwarfs issued from mergers that failed to exceed the Chandrasekhar 


3-1Gyr = 10° years 
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limit represent precious observational endpoints to calibrate the intricate nu- 
merical simulations of binary interaction. 

Another field to which white dwarfs have a lot of exciting contributions 
to make is that of exoplanets. White dwarfs polluted by debris from their 
own planetary systems represent our only detailed window on the compo- 
sition of rocky planetary material outside our solar system. Through a de- 
tailed modeling of the atmospheres of those polluted objects and with a good 
understanding of element transport mechanisms in the envelope, it is possi- 
ble to precisely trace back the elemental composition of the infalling rocky 
bodies. So far, most objects have been found to have bulk compositions 
similar to that of rocky bodies in the solar system. Among the most inter- 
esting cases discovered to date are systems with water-bearing planetesimals 
[84, 85, 86, 87, 88], differentiated bodies (i.e., bodies large enough to form a 
core, mantle, and crust) [89, 90, 91, 92, 93, 94, 95, 96, 97], and objects that 
have intriguingly high Li [98, 99, 100] and Be [101] abundances. 


3. Physical Processes in White Dwarfs 


Here we briefly describe the structure and cooling of a typical C-O core 
white dwarf. We refer to several reviews of white dwarf evolutionary models 
for a more detailed discussion |[102, 26, 103, 27, 2]. For the purpose of clarity, 
we define the nomenclature that we will use hereafter to refer to certain types 
of white dwarfs and models. Our choice is intentionally simplified from the 
actual practice but is adequate for our discussion. We divide the star in three 
regions: the core, which is made of C-O and constitutes 99% of the mass; the 
envelope, which includes the thin layers (the remaining 1%) of helium and 
hydrogen wrapped over the core; and the atmosphere which refers to the very 
thin outer layer that form the visible surface of the star. The spectral type, 
e.g. DA, DB, etc., is used in the strict sense of its empirical definition (Table 
1) and does not imply a particular atmospheric composition. Because the 
spectral type is not a proxy for the composition of the atmosphere (although 
they correlate strongly), the latter will be described as H- or He-atmosphere 
when it is dominated (including pure composition) by H and He, respectively. 
Finally, to distinguish white dwarfs that have only an outer He layer from 
those with both a He and H layers in their envelope, we will use H-rich and 
He-rich white dwarfs, respectively, even though 99% of the star is made of 
C-O. Note that the composition of the atmosphere is usually but not always 
the same as that of the outermost layer of the envelope. While these terms 
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have strong overlap among white dwarfs, they are not interchangeable, hence 
the need for these distinctions. 


3.1. Evolution 


The evolution phases that follow the main sequence (red giant, horizontal 
branch, AGB, and post-AGB) represent a very small fraction of the lifetime of 
a star [104, 105], but they set most of the initial conditions for the subsequent 
white dwarf evolution, most importantly the mass fractions of H, He, C, and 
O and of traces of heavier elements, and the relation between the mass of the 
newborn star and its final mass as a white dwarf* [106, 107, 108, 109, 110, 
111]. These properties depend themselves on stellar rotation, the efficiency 
of convective mixing, nuclear reaction rates, and the metallicity at the time 
of star formation [112, 113, 109, 111]. 

While it takes only ® 3.5Myr for a white dwarf to cool from Tig ~ 
120000 K to 40000K [2] (compared to a full lifetime of several Gyr) this 
early phase of the evolution has been studied extensively due to the ease of 
observing these bright, hot, young stellar remnants. The high luminosity and 
negligible nuclear energy production imply very fast cooling rates. Residual 
nuclear burning is negligible in hot white dwarfs formed from solar-metallicity 
progenitors, but is predicted to be a moderate source of energy at very low 
metallicities [114]. 

Element sedimentation works very fast in young white dwarfs, resulting 
in a chemical stratification of the envelope and nearly pure-H atmospheres 
in less than 1 Myr. In contrast, H-deficient stars have atmospheres showing 
a mixture dominated by He, C, and O in roughly comparable proportions 
for up to 3 Myr. This surface composition is determined by the balance 
between stellar winds and gravitational settling. The top He and H layers 
are initially non-degenerate and expanded by thermal pressure. As these 
layers cool from Tyg = 10° to 40000K, the white dwarf shrinks by about 
60% [2]. This gravitational contraction is a minor source of energy, given the 
small mass of the outer layers. 


4This is the “initial-to-final mass relation.” The final mass is always smaller due to 
extensive mass loss during the life of the star, particularly during the post-main sequence 
phases. 

5In astrophysical nomenclature, which can be arcane, “metals” refers to all the elements 
heavier than H and He. The metallicity of a star is its mass fraction of Z > 2 elements, 
which is ¥ 0.013 for the Sun. 
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Figure 4: Cooling track illustrating the behavior of the neutrino luminosity D, in relation 
to the photon luminosity L, in a representative DA model of 0.6 Mo. The star cools from 
left to right. For the first 20 Myr, the star loses more energy by neutrino emission than 
by radiation. This figure is based on the evolutionary tracks of Ref. [2]. 
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Neutrino cooling is an important energy sink during the early stage of 
evolution [27], where the interior is hot enough for neutrinos to be created by 
electroweak interactions. The vast majority of the neutrinos escape directly 
from the central regions to the stellar exterior without interacting with the 
stellar plasma, thus resulting in a net energy loss. Figure 4 demonstrates that 
neutrino processes dominate the cooling time scale until the surface reaches 
a temperature of © 25000 K corresponding to a cooling age of 20 Myr. 

The next phase of evolution is thermal cooling, where the structure de- 
pends on the properties of the degenerate electrons and thermal ions. The 
slow leakage of thermal energy from the ions in the plasma largely sets the 
cooling rates. The interior initially acts as a moderately coupled plasma, 
but Coulomb interactions grow as the star cools and the plasma eventually 
crystallizes through a first-order phase transition (Section 3.2). 

The He and H envelope contains relatively little mass, but plays an im- 
portant role in the cooling of a white dwarf because its large radiative opacity 
regulates the flow of energy to outer space. The final phase of cooling is when 
superficial convection, initially restricted to the photosphere, grows deeper to 
reach the degenerate interior, an event named convective coupling (Section 
Seay 

Figure 5 illustrates the thermal phase of the cooling process for evolu- 
tionary models in the mass range 0.4 — 1.3 Mo and Tyg range 50000K — 
2000 K. It highlights the three important events in the lifetime of a white 
dwarf: the end of the neutrino cooling phase and the beginning of thermal 
cooling (small filled circles), crystallization in the core (blue region), and con- 
vective coupling between the surface and the nearly isothermal core (small 
squares). The red isochrones show the strong dependence of the cooling rate 
on the mass and effective temperature of the white dwarf and how the mass 
determines the timing of the three main events in white dwarf cooling. A 
striking feature of the cooling tracks is that they are very nearly linear in this 
log-log figure and almost parallel. Since a star’s luminosity L, = 47 R?0T4,, 
where o is the Stefan-Boltzmann constant, this indicates that they cool at 
very nearly constant radius R, and that the mass-radius relation is essentially 
fixed after a few tens of Myr of cooling, with only a few percent deviation 
from that of a fully degenerate star. Notably, the most massive white dwarfs 
cool fastest even though they have a higher thermal energy content and a 
smaller surface to radiate heat. This comes from the early release of energy 
due to crystallization at high luminosities and the subsequent Debye cooling 
(Section 3.2). 


It 
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Figure 5: Evolutionary tracks (thin grey lines) of C-O core DA white dwarf models (M, = 
0.4 to 1.3 Mo in steps of 0.1 Me). The red curves are isochrones, labelled with the cooling 
time in units of Gyr. The small filled circles indicate the transition between the neutrino 
cooling phase, which dominates at higher luminosities, and the thermal cooling phase. 
The blue section of each evolutionary track corresponds to the crystallization of the C-O 
core. The set of square symbols indicates the onset of convective coupling between the 
surface and the degenerate core. The stars evolve from left to right, hence the Tp axis is 
backwards. This figure is based on the evolutionary tracks of Ref. [2]. 
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In summary, the computation of cooling ages suffers from three main 
uncertainties. First, the constitutive physics (equations of state, radiative 
and conductive opacities, neutrino production rates, diffusion coefficients) 
impacts the three important events described above. The second uncertainty 
is the amount of ion thermal energy that is stored in the interior of the star, 
which depends on the total number of C and O ions. Their relative proportion 
follows from the nuclear fusion rates and convective mixing during the pre- 
white dwarf evolution. Finally, cooling rates also depend on how fast this 
energy is transferred from the interior to the surface across the thin but 
opaque envelope, which depends on the total residual masses of H and He, 
that are set by the post-AGB evolution. 


3.2. Crystallization 


The onset of crystallization, the phase transition from liquid to solid of the 
C-O plasma at the center, is indicated by the upper end of the blue sections of 
the tracks in Figure 5. Because the internal energy is discontinuous between 
the liquid and solid phases, this predicted phase transition is of the first order. 
This corresponds closely to the first order liquid-solid phase transition of the 
one-component plasma and releases a substantial + 3 kpT per ion of latent 
heat [115, 116, 117]. We note that because massive white dwarfs have larger 
internal densities (for comparable temperatures), they develop a crystallized 
core at higher effective temperatures and earlier in their evolution. It takes 
of the order of 1-3 Gyr for crystallization to grow and encompass >90% of 
the degenerate core (Figure 5). The amount of energy released during this 
process, including an associated increase in gravitational binding energy (of 
comparable magnitude to the energy released by latent heat, see Section 5.5), 
introduces a cooling delay of about 1-2 Gyr. 

It was demonstrated fifty years ago that the slow down of the cooling 
process due to the phase transition would result in an observable statistical 
excess of white dwarfs at certain temperatures and masses, corresponding 
to core conditions where crystallization occurs [115]. The first evidence of 
such an accumulation of white dwarfs was found in the cooling sequence of 
the globular star cluster NGC 6397, which was interpreted as the release of 
latent heat from crystallization [118, 119]. More recently, the very precise 
data from the Gaia space observatory revealed clear evidence of the early 
crystallization among white dwarfs with a mass higher than the average in 
the disk of the Milky Way (Figure 5) [120, 13]. 
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Theoretical studies of crystallization following the initial prediction char- 
acterised the crystallization process and associated physics. In particular, 
crystallization can induce phase separation, gravitational sedimentation and 
distillation in a mixture of C, O, Ne and traces of higher mass elements 
[121, 122, 123, 124]. The Gaia observations spurred new, detailed studies 
[14, 125, 126, 17, 127]. The physical processes associated with crystallization 
are linked to the atomic weight of several key elements and isotopes crystalliz- 
ing at different temperatures and inducing a change in chemical stratification 
that releases additional gravitational energy. These new developments are 
discussed in Section 5.5. 

Figure 5 shows that more massive models have a rapid final phase of cool- 
ing following core crystallization. This produces the extensive deformation of 
the later isochrones in the figure. This phase of “Debye cooling” occurs when 
the core temperature falls below the Debye temperature of the Coulomb solid 
6, = hwy, where w, is the plasma frequency of the ions in the lattice. At 
such low temperatures, the discrete nature of the phonon spectrum of the 
lattice becomes apparent and the specific heat of the solid, initially given 
by the classical high-temperature Dulong-Petit value of Cy = 3k, per ion 
now decreases as ~ (T'/05)°. In other words, the ions in the solid have lost 
most of their capacity to store thermal energy, which initiates this final and 
relatively rapid cooling phase [115, 24]. In white dwarfs with a typical mass 
of 0.6 Ms, Debye cooling starts very late in the evolution and the oldest such 
stars in our Galaxy have just recently entered the phase of Debye cooling. 


8.3. Energy transport 


There are three main mechanisms of energy transport in white dwarfs: 
radiation, conduction and convection. These methods of transport can occur 
in parallel. This is augmented by neutrino cooling in early stages, which 
can transport energy directly from the core to the outside. Conduction from 
degenerate electrons dominates energy transfer within the degenerate core 
throughout the cooling process, resulting in fairly isothermal cores. Convec- 
tion and radiation dominate in the non-degenerate portion of the envelope, 
which we now discuss. 

As they cool, models of H-rich white dwarfs with a thick H layer (My/M, ~ 
10~*) start to develop convective instabilities at the photosphere due to the 
partial ionization of H at an effective temperature of around 18000 K [129], or 
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log q = log|l — M(r)/M,] 
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Figure 6: Evolution of the convection zone and radiative opacity (blue contour surfaces, 
labeled with log Kk, in cm? g~') in the outer envelope of a 0.6 Mo H-rich white dwarf as 
a function of the fractional mass above a given layer in the star (a vertical coordinate) 
and Teg. The star cools from left to right. The convection zone is bounded by the black 
solid contours and is largely determined by regions of high radiative opacity that occur 
when H is partially ionized. The depth of the photosphere (defined as the level where 
the Rosseland mean optical depth 7, = 2/3; dashed line) and of the onset of electron 
degeneracy (0 = T/T; = 1, where JT; is the Fermi temperature; dotted line) are also 
shown. This figure was generated using outputs from the STELUM white dwarf evolution 
code [128]. 
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a cooling age of 100 Myr.°® Initially convection is inefficient and characterised 
by low Péclet numbers (see Figure 16 of Ref. [129]). In Figure 6, the con- 
vection zone eventually grows deeper to sub-photospheric layers at around 
12000 K where it becomes efficient at transporting energy (Péclet number 
of order unity). At Tog < 10000 K, the convective photosphere becomes in- 
creasingly more adiabatic’ and less dependent on the assumed efficiency in 
the 1D model of convection [130], while the convection zone continues to grow 
deeper where the dense flows are almost fully adiabatic. At the same time, 
the degeneracy front moves upward as the white dwarf cools (Figure 6). This 
creates the conditions for the so-called convective coupling event [102, 26] as 
the convection zone reaches the degenerate layers at a Tog of about 6000 kK, 
following a rapid change in envelope opacity. This is an important event in 
the evolution of a white dwarf since this allows, for the first time, a direct 
coupling between the energy reservoir in the degenerate interior with the 
radiating surface. The transport of energy from the degenerate core to the 
photosphere no longer depends on the opacity of an intermediate radiative 
layer. The efficient convection allows internal heat to escape rapidly, tem- 
porarily slowing down the decrease in luminosity. Eventually, the main effect 
is however to accelerate cooling. After convective coupling, the convective 
zone stops growing inward (Figure 6). Convective layers never reach deeper 
than the superficial My/M, ~ 10~° mass fraction. Compared to main se- 
quence stars, these relatively shallow convection zones never reach densities 
large enough to produce observable magnetic dynamo effects [131]. Hence, 
and unlike the solar magnetic field, the population of magnetic white dwarfs 
cannot be explained by a dynamo induced by an outer convection zone (Sec- 
tion 5.1.2). 

He-atmosphere white dwarfs follow a similar pattern of development. of 
convective instabilities due to the partial ionization of He. Convection ap- 
pears near the photosphere at ~ 50000K, the convection zone deepens 
rapidly at © 30000 K, and convective coupling occurs at © 12000 K [{102, 8]. 
Compared to H-dominated atmospheres, these processes are shifted to much 
larger effective temperatures owing to the higher ionization potential of He. 


Tn Figure 6, which is based on a simple 1D model of convection, the onset of convection 
occurs at 16500 K. 

?Convection is said to be adiabatic when a rising fluid eddy does not exchange energy 
with its surroundings before dissolving higher up in the star. In this limit, the vertical 
temperature-pressure profile corresponds to an adiabat in the equation of state. 
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Even though convection involves complex three-dimensional flow pat- 
terns, most stellar models rely on the mixing-length theory of convection 
(MLT), a simple 1D phenomenological model of the convective flux [132, 
104, 133, 105]. The MLT considers a fluid element that, due to some local 
perturbation, is displaced vertically. As it rises and moves down the star’s 
temperature and pressure gradients, it expands and cools, and suffer radia- 
tive losses to the surrounding fluid. If after this initial perturbation the fluid 
element is buoyant, it will continue to rise and the region is convectively un- 
stable. The MLT provides estimates of the convective velocity, the convective 
energy flux and the resulting temperature gradient in the convective region, 
all functions of a single free parameter determining the vertical length scale 
of convective flows. This parameter is known as the mixing length, which is 
typically 1-2 times the local pressure scale height Hp in the star, which is 
defined by 

P(r) « exp(—r/Hp) , (2) 


where, for an ideal gas in hydrostatic equilibrium, H, = P/(dP/dr) = 
kT /gu, where g is the gravitational acceleration and p the average mass 
per particle. Except in regions of low optical depth (r < 1), the fluid el- 
ement rises and expands nearly adiabatically. Such adiabatic convection is 
insensitive to the choice of mixing length and can carry very large fluxes. 
On the other hand, white dwarf model spectra are very sensitive to the ef- 
ficiency of convection at the surface, which is largely unconstrained by the 
MLT [134, 135]. Nevertheless, the MLT provides a good description of the 
importance of convection in the evolution of a white dwarf. Modern 3D 
radiation-hydrodynamics simulations of convection in white dwarfs are dis- 
cussed in Section 5.1. 

Convection is thought to play a major role in spectral evolution, as ver- 
tical mixing can change the composition of the atmosphere (Section 2.3). 
Here we focus on major transitions from H- to He-dominated atmospheres. 
Convective dilution is a process in which the underlying He layer becomes 
convective in white dwarfs with very thin outer H layers (My/M, ~ 107!°). 
This results in a full mixing of the H layer into the larger underlying He layer 
at Teg of 30 000-20 000 K, resulting in a transition from H-dominated to He- 
dominated atmospheres [136]. Another important process is convective mix- 
ing, where models indicate that the top convective H layer, if My/M, < 10~°, 
reaches the underlying He layer (Figure 6). This also results in the full mix- 
ing of the H layer with the much larger underlying He layer [8, 32, 33]. For 


23 


thicker H layers, the superficial H convection zone never reaches the under- 
lying He layer and convective mixing cannot happen at any point during 
the evolution. We note that the accretion of solid planetary debris can also 
induce more moderate changes of atmospheric composition, such as varia- 
tions in the H/He ratio in He-dominated atmospheres from the accretion of 
water-rich asteroids as a source of H [87]. 

The internal reservoir of thermal energy in a cooling white dwarf is in- 
sulated by the radiative H and He envelope. It is the vertically integrated 
opacity of this non-degenerate envelope that sets the rate of energy loss. Fig- 
ure 6 illustrates this opacity barrier for a DA white dwarf, highlighting the 
contribution of the insulating layers with large opacity (colored contours) 
above the degeneracy boundary which is defined as 8 = T/T, = 1, where 
Ty is the Fermi energy (dotted line). The opacity maximum is caused by 
the partial ionization of H. The upper boundary of the core’s reservoir of 
ion thermal energy roughly corresponds to the degenerate boundary, where 
the running radial integral of the specific heat [ CyT' dm reaches a plateau 
[137]. The strength of the insulating radiative layer is proportional to the 
optical depth t; = — { Krpdr between the base of the convection zone and 
the degeneracy boundary 0 = 1. Because density increases with depth, the 
integrand is sharply peaked in the layers just above the degenerate boundary. 
This demonstrates how this region controls the rate of energy transfer from 
the degenerate core to the exterior. The low-density atmosphere (tT ~ 2/3) 
provides a negligible contribution to the opacity barrier and does not in- 
fluence the cooling process. This implies that the cooling rates, hence the 
relation between core and surface temperatures, are exactly the same for H- 
rich and He-rich white dwarfs until the convection zone ultimately reaches 
the degenerate boundary. This convective coupling results in the sudden 
disappearance of this insulating layer. From now on, the degenerate core 
is connected to the surface by a nearly adiabatic convection zone that can 
carry far more flux to the surface. The opacity of the atmospheric layers 
now controls the cooling process, leading to age differences of & 1 Gyr be- 
tween H-rich and He-rich white dwarfs once they reach surface temperatures 
of 4000 K [2]. 

Conductive opacities in the deep interior are relatively well established 
[138]. A better accounting of electron-electron scattering has resulted in an 
improved understanding of these opacities in the regime that corresponds to 
the core-envelope interface that largely governs cooling rates [139]. Conduc- 
tion works in parallel with radiative transfer at the base of the insulating 
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layer. In optically thick layers, radiation transport is described as a diffusive 
process and only depends on the Rosseland mean opacity, which is rather 
well known in this regime (see Section 4). 

The atmosphere does not impact cooling rates for most white dwarfs, but 
important work has also been performed to understand radiative transfer, 
its effect on line formation, its interaction with convection, as well as the 
radiative levitation of trace metals. In particular, non-local thermal equilib- 
rium effects have a major impact on the structure of the upper atmosphere 
in young hot white dwarfs (Tog 2 40000 K) [140, 141]. However, apart from 
the upper atmosphere of hot white dwarfs, local thermal equilibrium prevails 
in the entire star. 


3.4. Diffusion 

Atomic diffusion is a key process in our understanding of white dwarf 
structure and evolution. It is responsible for the relaxation of the com- 
position profile resulting from pre-white dwarf evolution into a chemically 
stratified young white dwarf with a C-O core surrounded by He and H lay- 
ers (Figure 1). In competition with stellar winds and radiative levitation, 
diffusion plays a central role in the subsequent spectral evolution [2]. Even 
once the white dwarf has cooled below 30000 K and the internal structure 
has largely settled, atomic diffusion couples to new mechanisms that appear 
later in the evolution, such as convection, accretion and crystallization. In 
DZ white dwarfs, it differentially sorts the infalling planetary material and 
alters the observed elemental composition at the photosphere. The upward 
tail of carbon at the core/He transition (Figure 1) can be dredged up by 
the He convection zone to “pollute” the surface of DQ white dwarfs with 
traces of carbon. During crystallization, the associated fractionation of trace 
elements between the solid and liquid phases couples the phase transition to 
the release of gravitational energy and diffusion in the liquid phase (Section 
5.5.2). 

Diffusion in stellar plasmas occur in the presence of a gravitational field 
and a radial temperature gradient. The inter-diffusion of ions is governed 
by the combined effects of gradients of concentration, pressure, and tem- 
perature, and an induced electric field. Composition transition zones are 
the regions of continuously varying composition that separate homogeneous 
layers. For instance, a white dwarf with a H layer has two composition tran- 
sition zones, between the H and He layers, and between the He layer and 
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the C-O core. Although this is not necessary, most evolution models as- 
sume that the abundance profiles in the transition zones are specified by the 
condition of diffusive equilibrium and neglect thermal diffusion. The result 
is interpenetrating diffusion tails at the transition zones, as can be seen in 
Figure 1. 

Other phenomena in white dwarfs arise from the diffusion of a trace ele- 
ment in a (usually) fully ionized background of H and He. In the presence of 
surface convection that maintains a homogeneously mixed composition, dif- 
fusion allows heavier elements to drop out at the bottom of the mixed region. 
The diffusion of trace elements usually requires time-dependent calculations 
(6, 142, 143, 9, 144]. These calculations show that the e-folding diffusion time 
scale at the bottom of the convection zone ranges from days to ~ 10‘ years 
[145, 144] and that diffusive equilibrium is reached in a matter of at most 
millions of years in white dwarf envelopes but can take much longer in the 
degenerate interior [146, 136]. The short diffusion time scale near the surface 
is an important result as it implies that metal-polluted DZ white dwarf at- 
mospheres must have accreted these elements very recently in their evolution 
or are currently accreting, as their sinking times are short compared to evo- 
lutionary time scales. Otherwise, they would have long disappeared below 
the photosphere. This led to the current understanding where these elements 
result from recent or continuous accretion of planetary material—perhaps in 
the form of asteroids that are tidally disrupted—onto the white dwarf where 
they can be observed in the spectrum [39]. Conversely, the upward diffusion 
of small amounts of H in the He layer can counteract convective dilution, 
where the He convective layer reaches up to the H/He transition zone and 
otherwise results in a mixed composition layer [136]. 


3.5. Pulsations 

The properties of the H-atmosphere (DAV) and He-atmosphere (DBV) 
pulsating white dwarfs are closely linked to the appearance of convection, 
partial ionization zones and large peaks in opacities. The opacity bumps due 
to the photoionization of H or He become so large that convection sets in 
and the pulsation modes are driven (linearly unstable) mostly through the 
convective driving mechanism. This comes from the fact that the integrand 
dW/dm of the work integral® [147, 148], which gives the contribution of 


’The kinetic energy of a pulsation mode is sustained by the work done by both the 
pressure and the gravitational forces averaged over a pulsation cycle. 
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each stellar layer to the net driving/damping of a mode averaged over one 
oscillation period, has most of its driving contribution near the base of the 
superficial convection zone [56, 58]. Because the convective turnover time 
scale at the base is much shorter than the pulsation periods of interest, the 
convective flux responds almost instantaneously to pulsations [149]. 

The compact nature of a white dwarf implies that the period of a pulsation 
mode with eigenfunction of radial order k and degree index @ is much smaller 
than in other stars. The observed modes are all internal, non-radial gravity 
waves, with periods © 100-2500 seconds. The eigenfunctions of these g modes 
have large amplitudes mostly restricted to the external envelope layers, the 
exact opposite of what is seen in non-degenerate stars like the Sun where g 
modes probe the deep interior. 

The condition for local g mode propagation is that the square of the 
frequency of a mode has to be less, at a given depth in a model, than the 
smaller of the square of the Brunt-Vaisala frequency 


CydinP dilnp 
N? = = 
" (G dr dr 8) 
and the Lamb frequency 
(e+ 1c 
gp = Ee (4) 


where c, is the sound speed, and C\ and C> are the specific heats at con- 
stant volume and pressure, respectively. The Brunt—Vaisala frequency is the 
oscillation frequency of a perturbed fluid element in a convectively stable 
region and the Lamb frequency is the inverse of the sound crossing time over 
the horizontal pulsation wavelength. The large values of these characteristic 
frequencies in white dwarfs, allowing for propagation of gravity modes in the 
envelope, is a direct consequence of electron degeneracy on the Brunt—Vaisala 
frequency [56]. 

The fact that g modes show maximum amplitudes in the outer non- 
degenerate envelope implies that the observed pulsations are most sensitive 
to the parameters of these layers, such as the extent of the convection zone 
and of the H and He layers. The composition transition zones are of crit- 
ical importance for asteroseismology. Their influence is revealed by peaks 
in the Brunt—Vaisala frequency at composition transitions that “pinch” the 
eigenfunction of the pulsation modes, leading to spatially confined “trapped” 
modes and resonances. Diffusion determines the sharpness of the transition 
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zones and their potential for mode trapping. Because of this phenomenon, 
asteroseismology provides a powerful method to constrain the internal strat- 
ification of white dwarfs and to determine the thicknesses of their He and H 
layers [56]. In contrast, the degenerate interior of a white dwarf is difficult 
(but not impossible [73, 74]) to survey using pulsations since the eigenfunc- 
tions have very low amplitudes in those regions. 

In recent years, progress has been made to understand the boundaries of 
the DBV and DAV instability strips in terms of Tyg and surface gravity. A 
current problem in white dwarf asteroseismology is that the predicted Tog of 
the lower boundaries are lower than the observed values by a factor of about 
two [150, 151]. It has recently been discovered that DAV pulsators near the 
observed low-Tug edge ( 10500 K) exhibit energetic outbursts that increase 
the luminosity by ~ 10% on a non-periodic time scale of days [152, 153]. The 
current explanation for this behaviour involves a resonant transfer of energy 
between driven and damped modes [154]. 


3.6. Magnetic fields 


The origin of magnetic fields in single white dwarfs is still largely unex- 
plained and we refer the reader to several reviews on the topic [22, 155, 156]. 
Most of the major scenarios involve external processes or pre-white dwarf 
evolution rather than an origin during white dwarf cooling. They include 
binarity and common-envelope evolution, post-main sequence or main se- 
quence internal dynamos, flux conservation of fossil fields from the pre-main 
sequence, and double white dwarf mergers. One additional proposed sce- 
nario, consisting of an internal convective dynamo generated during the 
crystallization process, could happen naturally during the cooling process 
[157, 158, 159]. However it is still unclear if this scenario, whose efficiency 
depends on the rotation velocity and the diffusion of internal magnetic fields 
to the surface, is a significant channel for the creation of magnetic white 
dwarfs [156, 160]. 

There is evidence that magnetic field incidence increases with mass [22, 
21, 156] —a trend that can be explained by all scenarios for the origin of 
magnetic white dwarfs. Recent studies have shown that all magnetic field 
strengths in the range 0.1-500 MG are equally likely [156]. There is also 
evidence that younger white dwarfs (<0.5 Gyr) have a lower incidence of 
magnetic fields and smaller field strengths [156, 160]. Some of these changes 
with cooling age could be caused by the outward diffusion to the surface 
of a magnetic field generated during the pre-white dwarf phases or from a 
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crystallization-induced dynamo, which must overcome Ohmic decay, with 

estimated time scales of 0.5-4 Gyr depending on white dwarf mass [131]. 
Magnetic fields can affect the structure (e.g. the mass-radius relation) of 

white dwarfs if the magnetic pressure is larger than the central pressure, 


= > P.~ 104 erg em? (5) 
87 : e 

or B => 5 x 10"G. While this simple estimate only accounts for the direct 
contribution of the magnetic field to the total pressure and ignores its effect 
on the equation of state (EOS), it is consistent with more detailed calculations 
that show that the mass-radius relation is not affected unless the interior field 
reaches ~ 10'° — 10'*G [161]. Such fields are orders of magnitude higher 
than the strongest fields observed at the surface and are considered unlikely 
in single white dwarfs (Section 4.4). 

Similarly, observed magnetic fields are too weak to influence white dwarf 
cooling. This is because the bottleneck in cooling rates is at the interface 
between the degenerate core and the envelope above. Magnetic fields smaller 
than 1000 MG are not strong enough to impact energy transfer rates in 
these deep layers. The interaction between convection and magnetic fields is 
discussed in Section 5.1.2. 

Magnetic fields in white dwarfs are detected and characterized by their 
effect on their spectra through Zeeman splitting and shifts of spectral lines 
or polarization of the continuum. These effects, which can be quite strong, 
are surface manifestations of magnetic fields where the density and temper- 
atures are low compared to the interior. The problem of modelling spectral 
line profiles in the presence of comparable Stark and Zeeman effects has 
remained largely unsolved until very recently [162]. This means that spec- 
tral lines of strongly magnetic white dwarfs cannot yet be used for reliable 
atmospheric parameter determinations or for studying magnetic effects on 
the atmospheric structure. Work is ongoing for calculating line splitting for 
highly magnetic white dwarfs with field strengths above the Paschen—Back 
regime, with the goal of better understanding magnetic field geometries and 
atmospheric composition [163]. On the other hand, analysis with the pho- 
tometric method using the accurate distances obtained with the Gaia space 
observatory is a reliable alternative to determine the stellar radius of the 
cooler magnetic white dwarfs [21] rather than fitting the gravity from the 
observed line profiles (Section 2.3). 
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3.7. General relativity in white dwarfs 
The importance of general relativistic effects on the structure of white 
dwarfs can be estimated from the ratio 
R, 2GM, 
= — = 6 
a Ae c?R, ( ) 
where R, is the Schwarzschild radius. For a typical white dwarf of 0.6 Mo 
with a radius of 8900 km, y = 2x 10~4 and the effect is commonly neglected in 
white dwarf models. However, a massive white dwarf near the Chandrasekhar 
mass limit of 1.4 Mo with a radius of & 1000km has y = 0.004. Thus, the 
most massive white dwarfs are compact enough for weak general relativistic 


effects to affect their structure. In particular, this reduces the Chandrasekhar 
mass by 2-3% [164]. 


4. Physical regimes in white dwarfs 


In this section, we examine the physical regimes encountered in white 
dwarfs as a reference for the development of models, theories and experiments 
to address the various challenges they pose. The physical conditions inside 
white dwarfs vary over many orders of magnitude from the very dense, fully 
ionized C-O core to the tenuous gas typical of the visible photosphere. These 
conditions also evolve as the star cools. It follows that a wide range of 
physical regimes are represented. 


4.1. Models with an outer hydrogen layer 

Figure 7 shows a series of density—-temperature profiles of a representative 
model for a white dwarf with an outer H layer with total mass M, = 0.6 Mo 
(see Figure 2) as it cools from Tyg = 30000 to 3000K.° Such white dwarfs 
are of spectral types DA and turn into spectral type DC below Tyg 5000 K 
(Section 2.3). The model assumes a core with equal mass fractions of C and 
O surrounded by a He layer of 107? M, and an outer H layer of 107-4 M,. As 
is justified for the vast majority of single white dwarfs, rotation and magnetic 
fields are neglected. The 30000K model has a white dwarf cooling age of 
about 9Myr and the coolest and the oldest white dwarfs known have Tig 
somewhat below 4000 K and an age of about 10Gyr. This set of models 
spans all but the very earliest cooling history of H-rich white dwarfs. 


°The models presented in Sections 4.1 to 4.3 are computed with the white dwarf evo- 
lution code STELUM [128]. 
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Figure 7: Density-temperature profiles of a cooling H-rich white dwarf model. Profiles 
from Tyg = 30000K to 3000 K in steps of 1000 K are shown. The cooling sequence runs 
from top to bottom. The thickness of the line indicates the dominant atomic constituent, 
and regions that are convectively unstable are highlighted in green. The different shadings 
of blue indicate where the gas is mostly composed of Hz, atomic H gas, or a partially 
ionized H plasma. The white background indicates a fully ionized plasma of C, O, He 
or H. Finally, the circles on each profile mark the location of the photosphere, and the 
dotted lines indicate the values of some important plasma parameters (see text). Those 
profiles were obtained assuming a core with equal mass of C and O, an envelope with 
log(Mye/M,) = —2 and log(My/M,) = —4, and a total mass of M, = 0.6 Mo. In this 
log-log plot, the outer layers of the star are greatly expanded. The surface of the star is 
toward the left, the center is on the right. 
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The stratified structure of the star in C-O, He and H layers is revealed by 
the thickness of the grey lines in Figure 7. The C-O core is nearly isothermal 
due to the high thermal conductivity of the degenerate electrons. As the 
star cools, the central temperature decreases from 6 keV to 40 eV, but since 
it contracts very little during this part of the evolution, the density profile 
of the C-O core remains unaffected, ranging from log p(gcm~?) ~ 4.6 to 
6.6. Under these conditions, both C and O are fully ionized. The He layer 
is also fully ionized and mostly isothermal. The temperature and density 
decrease by several orders of magnitude in the thin H layer, all the way 
to the photosphere (solid dot). Initially fully ionized in the hotter models, 
the H layer enters a regime of partial ionization around Tyg = 15000K. 
At the photosphere, H atoms dominate below 10000K and H»2 molecules 
form in models cooler than 5000 K. Pressure dissociation and ionization of H 
occurs in the coolest models at p ~ 1gcm7?. Energy transport towards the 
surface is dominated by electron thermal conduction and radiation, but for 
Tog S 18000K a convectively unstable region develops and grows in the H 
layer (thick green line; compare to Figure 6). 

The (p,7) plane can be divided into regions of various plasma regimes. 
Electron degeneracy sets in when the degeneracy parameter, the ratio of the 
temperature to the Fermi temperature given by 


kel 2kaIm, 


is smaller than unity. Here, k, is Boltzmann’s constant, €, the Fermi energy, 
m, the mass of the electron, h the Planck constant and n, is the electron 
number density. Two lines of 0 = 1 are shown in Figure 7 corresponding 
to pure H and pure He compositions. The corresponding line for the C-O 
plasma is the same as that for He. Electrons are degenerate in the entire C-O 
core, the He layer and in the deeper regions of the H layer. At the observable 
photosphere, electrons are always classical and recombine to form atoms and 
molecules in the cooler models. 
The ion plasma coupling parameter 


Loe 
— akpT (8) 


is the ratio of the electrostatic potential energy of two ions of charge Ze at a 
separation of an ion sphere radius a over their kinetic energy (~ kgT). Ions 
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are strongly coupled and non-ideal when I’ > 1. Lines where [ = 1 are shown 
for pure H and pure He compositions. The corresponding line for the C-O 
plasma lies well above the [,. = 1 line and is omitted for clarity. The H layer 
is usually weakly coupled except in models cooler than 5000 K. Similarly, the 
He layer is initially weakly coupled for models hotter than + 20000 K but it 
eventually becomes strongly coupled. The plasma coupling is always strong 
throughout the C-O core. It is well known that a strongly coupled plasma, 
more specifically the one-component plasma model, freezes at [ & 175 [165]. 
This freezing line for the C-O plasma intersects the core structures for models 
with Tx¢ < 6000 K. The entire C-O core has become solid once the star has 
cooled to ¥ 4000 K (Section 3.2). 

At the extreme densities found at the center of white dwarfs, the Fermi 
energy of the electrons becomes comparable to their rest mass energy. The 
relativity parameter is defined as 

— Pr _ A(3n2n,)1/8 . (9) 


MeC Mel 


where py is the Fermi momentum and c is the speed of light. The z,,¢o = 1 
line indicates that electrons are moderately relativistic at the center of white 
dwarfs throughout their evolution. 

Finally, ions can display quantum mechanical behavior at very high den- 
sities and low temperatures. The !?C and !°O nuclei are bosons of spin zero 
and are subject to Bose-Einstein statistics. Quantum mechanical effects be- 
come manifest when the thermal de Broglie wavelength A, of a particle is no 
longer negligible compared to the interparticle separation. In this case, ions 
no longer interact as classical point particles. This effect is measured by the 


quantum parameter 
1/2 
Qrh? 


Wikked: 


The line of Ago = 1 is below all the white dwarf structures shown in Figure 7 
but we can expect that weak quantum effects occur in the solid C-O core of 
the coolest white dwarfs. 

Figures 8 — 10 offer a closer look at key plasma parameters and prop- 
erties throughout the models. As in Figure 1, they use a logarithmic mass 
coordinate log q for the radial position in the star 


log q = log (: _ “0 ' (11) 


Mth 
a a 


(10) 


* 
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where m(r) is the mass enclosed in a volume of radius r. This mass coordinate 
greatly expands the outer part of the star where much interesting physics 
occurs. As a point of reference, the photosphere (“surface”) of the star is 
located around logqg ~ —17 for the hotter models and can go as deep as 
log gq = —12 at Tig = 3000 K (Figure 6). 

The mass fractions of C (black), O (red), He (blue) and H (magenta) are 
shown in Figure 8a. As a result of the processes that occur in pre-white dwarf 
evolution and subsequent gravitational settling (Section 2.1), the structure 
of the star is stratified in well-defined layers according to composition. The 
composition boundaries are not sharp however, which reflects the diffusive 
equilibrium between the species (Section 3.4). In the subsequent panels of 
Figure 8, the regions where conductive or radiative transport dominate are 
shown in blue, the convection zone is in red and the central crystallized region 
in green. A quantity that is useful in interpreting subsequent figures is the 
average ion charge throughout the star (Figure 8b). In the fully ionized C-O 
core, the average charge (Z) = 6.86. It drops to (Z) = 2 in the He layer 
and then to (Z) = 1 in the outer H layer. Throughout most of the star 
(logq > —8), (Z) is unaffected by the star’s cooling. In the outer reaches 
of the cooler models, H recombines and the gas becomes neutral ((Z) — 0). 
Convective energy transport (in red) is clearly associated with the partial 
ionization of H. 

The pressure profiles of all the models in this sequence are nearly identical 
(Figure 8c). This is because the EOS is dominated by the degenerate electron 
pressure with O < 1 and is insensitive to the temperature in the cooling star. 
Figure 8d shows a drop in density at the He/H boundary (log g = —4) but 
there is no corresponding feature at the boundary between the C-O core and 
the He layer (log g = —2). The requirement that the pressure be continuous 
in a star in hydrostatic equilibrium implies that P. (* P) and thus n, all be 
continuous in a highly degenerate plasma. This feature in the mass density 
profile appears when Z/A changes from 0.5 to 1 at the He/H boundary. 
Figure 9a shows the same temperature profiles as Figure 7 but in terms of a 
different radial coordinate. 

The ion plasma coupling parameter (Equation 8) is quite large in the 
C-O core, with values of [ = 50 in the hotter models to above 1000 in the 
crystallized core of the 3000K model (Figure 9b). The drop in (Z) from 
6.86 in the C-O core to 2 in the He layer reduces [ by roughly a factor 
of 10, but the He plasma remains strongly coupled except in the hotter 
models. The strength of the coupling decreases again upon entering the H 


34 


8 | aap ian Garam | ams aa Gama | a a ima & ] T 


i rs : L 1] 
[ 7 C b 4 
ra 0.8 & a” 6E 4 
ie} z El 
3 IE il L J 
3 0.6 = ae 1 
ce L J N 46> + 
E 7 Vv L 4 
0.4 + L 4 
o r a L I 
E E al ae al 
0.2F 4 L 1 
0 L | @) lem | 4 Ne | it | a i | 
0) 2 4 6 8 10 @) 5 10 15 
log (1-m(r)/M.) log (1-m(r)/M,) 
25 oe | i im a | T | Lia a | al | ae oo T T | ae T | ama ior ] T 
c 7 af d 
P 20 - = i i 
[ al ; L 
3 i | = L 
1S) 
20 L 4 L 
A, [ J Q [ 
| | Q0 L 
cate j £ f{ 
iB N -5F 
5 1 1 1 1 | i 1 i i | 1 L 1 1 | 1 1 1 1 1 | L L i 1 | 1 1 1 1 | 
0) —5 —10 -15 0 9) 10 15 
log (1—m(r)/M.) log (1—m(r)/M.,) 


Figure 8: Properties of a cooling sequence of models of a H-rich white dwarf of M, = 
0.6 Me consisting primarily of a C-O core, surrounded by a 107? M, He layer and an 
outer H layer of 1074 M,. Models range from Tyg = 30000K to 3000K (heavier line) in 
steps of 1000 K. The radial coordinate is a logarithmic mass coordinate (Equation 11) with 
the center of the star located on the left at logqg = 0 and the surface towards the right 
(note the “backward” X-axis). In panels b-d, conductive and radiative regions are in blue, 
convectively unstable regions in red and the crystallized core in green. a) Composition 
profile showing the mass fractions of C (black), O (red), He (blue) and H (magenta), b) 
Average ion charge (Z), c) Pressure (10/7 erg cm~? = 10/2 dyncm~? = 1 Mbar = 100 GPa), 
d) Density profile. 
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Figure 9: Same as Figure 8 but for the a) Temperature profile (leV = 11605K), b) Ion 
coupling parameter [ (Equation 8), c) Electron degeneracy parameter 0 (Equation 7), 
and d) Ion coupling and degeneracy parameters. 
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layer (logg < —4). The coupling is weak but not negligible in the H layer of 
most models, with [T = 0.05. The base of the H layer of the cooler models, 
including in the convection zone, is strongly coupled. The ion coupling drops 
rapidly once the plasma recombines to form H atoms. 

The most well-known characteristic of white dwarf stars, that they are 
degenerate stars, is illustrated in Figure 9c where the electron degeneracy 
parameter © (Equation 7) is less than unity in every model through nearly the 
entire star and can be as small as ~ 10~4 at the center. The degenerate region 
encompasses the entire C-O core (logg > —2), the entire He layer (—4 < 
log q < —2) and most of the mass of the H layer. A notable feature is that 
O is nearly constant in the convection zone (in red), owing to the identical 
T(p) dependence of the constant specific entropy profile in the adiabatic 
convection zone and that of a line of constant 0. Degeneracy is lifted (and 
O rises rapidly) when the H plasma recombines near the surface. 

From the point of view of plasma physics, it is interesting to look at 
white dwarf structures in terms of [ and © simultaneously (Figure 9d). The 
models fall along a broad diagonal band that spans from strongly degenerate 
and strongly coupled plasmas ([, © > 1, lower right quadrant) to the weakly 
degenerate, weakly coupled plasmas ([ < 1, © < 1, upper left quadrant). 
Plateaus of nearly constant O along each curve reflect the rapid change in 
(Z) at the composition boundaries (Figure 8a). In these models typical 
of the most common white dwarfs, the plasma does not enter a regime of 
strong coupling and weak degeneracy (upper right quadrant) but it crosses 
the regime of warm dense matter ([, © ~ 1) around the He/H transition 
region. 

Carbon and oxygen ions display weak quantum mechanical behavior in 
the core with a quantum parameter A = 0.02 that later rises steadily to 
moderate levels in the crystallized center of the coolest white dwarfs where 
A ~ 0.6 (Figure 10a). Jumps in A occur at the composition transitions 
between layers where the ion mass decreases and the thermal wavelength 
increase. Remarkably, ionic quantum effects are non-negligible (A > 0.01) in 
all models cooler than Tyg = 10000K everywhere except in the very outer 
regions. This includes all of the He layer, nearly all the mass in the superficial 
H layer, and most of the convective zone of these cool models. As a rule of 
thumb, a perturbative treatment of the ion quantum effects in the EOS, 
such as the Wigner—Kirkwood expansion to order h?, is usually adequate for 
A<0.1. 

Figure 10b shows the relativity parameter for the electrons (Equation 9). 
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Figure 10: Same as Figure 8 but showing the a) Ion quantum parameter A (Equation 10), 
b) Electron relativity parameter x, (Equation 9), c) Total opacity Kio, (Equation 15), and 
d) Adiabatic temperature gradient Vaaq (Equation 18). 
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Relativistic effects are important at the center of 0.6 Mo white dwarfs as x, ~ 
1 but decrease steadily outward. Relativistic corrections are not negligible 
(x, > 0.01) throughout the entire C-O core, the full He layer and most of the 
H layer (logg > —7). In more massive white dwarfs, the increased central 
density causes 2, to rise significantly (Section 4.3). 

Energy transport in stars can proceed through radiation, conduction or 
convection. Except near the photosphere, radiation transport can be de- 
scribed very accurately as a diffusive process of the radiation field in local 
thermodynamic equilibrium, i.e. a black body distribution at the local tem- 
perature T(r). The flux of radiation through a layer of radius r in the star 
can be written as [104] 


aT 1l6aT? dT 
Pes = — Dra, ied ) 12 
a(r) Var ohap ar ve 


where D,aq is the radiative diffusivity. This defines the frequency averaged 
Rosseland mean opacity Kp which is effectively an inverse diffusivity. Simi- 
larly, the conductive flux of is given by 
dT 
ser = = oon rae 13 
1 = —Deona (13) 
where Deon is the thermal conductivity. By analogy with Equation 12, we 
define a conductive opacity 


160T? 
DD eondp 


Keond = 


(14) 


As radiative and conductive transport can occur in parallel, the total opacity 
Ktot is 


1 1 1 
= ; (15) 
Ktot KR Kcond 
and the heat flux through a non-convective layer in the star becomes 
160T? dT 
FPr)y=- Sr 16 
age a aa (16) 


The high degree of electron degeneracy in the white dwarf core favors conduc- 
tive energy transport which eventually is superseded by radiative transport 
closer to the surface. The conductive opacity Kong is very small (very high 
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thermal conductivity) in the C-O core and in the He layer for the coolest mod- 
els (Figure 10c) but it increases rapidly towards the surface. The Rosseland 
mean opacity (Kg) is many orders of magnitude larger than the conductive 
opacity in the C-O core and is roughly constant outward until the H plasma 
starts to recombine and k, has a broad peak followed by a rapid decline near 
the surface where H is (nearly) all in atomic or molecular form (see Figure 
8b). The net result is that the total opacity Kio: (Figure 10c) is very small 
in the C-O core and He envelope, has a broad peak in the H layer for models 
with Teg < 12000 K, and drops rapidly near the surface. A striking feature 
of Figure 10c is that the region of convective energy transport (in red) coin- 
cides with the broad peak in k,,,. This is a consequence of the Schwarzschild 
criterion for convective instability in stars [105, 104] 


Vad < V ’ (17) 
where 
OlnT 
Vad = OlnP M8) 
iS 


is the adiabatic temperature gradient,!° a thermodynamic quantity obtained 
from the EOS, and 
OlnT 
~ OlnP 
is the dimensionless temperature gradient in the star. With some simple 
approximations that are justified in white dwarfs, the convective instability 
criterion can be rewritten as 


3 Ly, Pot P 
ec Kot 5 
64r0G M, T4 qe a2 


where G is the constant of gravitation and L, the luminosity of the star. 
Thus, aside from the factor P/T“, convective instability is favored when Vaa 
is small or Aitot is large. 

Figure 10d shows that the convective region generally occurs when V aq 
becomes small. Generally, in (p,7’) regimes where partial ionization occurs, 
Vaa decreases (see Figures 10 and 17 of Ref. [166] for example), Kona is very 
large because of the reduced electron degeneracy (hence, Kiot ~ Kp), and Kz 
increases due to the contribution of bound-free absorption to the opacity. 


(19) 


Vad < (20) 


10The derivative is taken at constant specific entropy S. 
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4.2. Models with an outer helium layer 

Internal structures of white dwarf models with only an outer He layer are 
shown in Figure 11. These models are representative of about 20% of all white 
dwarfs. The models shown have otherwise the same parameters as those in 
Figure 7. The physical properties of these He-atmosphere white dwarfs are 
very similar to those with a H atmosphere except near the surface where 
He has replaced H. A thin convection zone appears near the photosphere at 
Teg © 50000 K and becomes progressively deeper as the star cools. Pressure 
ionization of He (lighter blue shadings) occurs at densities of a few gcm7~? 
and T S 20 eV, conditions found in the convection zone of the cooler models 
shown and around the photosphere of models of 5000 K or less. Generally, 
electron degeneracy and ion Coulomb coupling become significant at higher 
Teg (earlier in the evolution) than in H-rich models. The photosphere of the 
coolest models is found at densities approaching 1 gcm~* where we expect 
strong non-ideal effects to leave an imprint on the spectra of those stars 
(Sections 4.5, 5.2 and 5.3.2). 

The simple composition profile of these models is shown in Figure 12a. 
The behavior of (Z) within the C-O core and the base of the He layer is 
the same as in the DA models (Figure 12b). The He layer remains fully 
ionized from logqg = —2 to —13 in the hotter models, where it then drops 
to a plateau where He remains singly ionized ((Z) ~ 1). In the cooler 
models, recombination of He starts steadily deeper in the model and decreases 
smoothly to neutral He. 

The pressure, density and temperature profiles are nearly identical as 
those of H-atmosphere models shown in Figures 8c,d and 9a. However, He 
has a lower opacity than H at low temperatures, resulting in a different 
atmospheric structure (Figure 16). The latter acts as the surface boundary 
condition for models of the interior and the relative transparency of the He 
atmosphere leads to faster cooling. The main effect is a slightly lower central 
temperature (less than a factor of 2) after convective coupling has occurred 
(Tee < 12000 K) and a larger crystallized region at a given Tog. 

The higher charge of fully ionized He leads to stronger ion plasma coupling 
and a higher degree of degeneracy in the upper reaches of the models than in 
the H layer of the H-rich models. These two plasma parameters are shown 
in Figure 12c. In addition to the shift to higher I in the upper reaches of 
the He layer (O = 1), a striking difference is that the convection zone (in 
red) of models with Tyg < 12000K extends well into the strongly coupled, 
strongly degenerate regime ([ > 1, 0 < 1, lower right quadrant). In all the 
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Figure 11: Same as Figure 7, but for He-rich white dwarfs. The blue shadings now indicate 
the ionization of He, see the legend. An envelope with log(My./M,) = —2 is assumed, 
and the white dwarf mass is 0.6 Mo. The surface of the star is toward the left, the center 
is on the right. 
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Figure 12: Properties of a cooling sequence of models of a He-rich white dwarf of 
M, = 0.6 Mo consisting of a C-O core, surrounded by a 107? M, He layer. Models range 
from Teg = 30000 K to 3000 K (heavier line) in steps of 1000 K. The radial coordinate is a 
logarithmic mass coordinate (Equation 11) with the center of the star located on the left 
at logq = 0 and the surface towards the right (note the “backward” X-axis). In panels 
b-d, conductive and radiative regions are in blue, convectively unstable regions in red and 
the crystallized core in green. Compare with the same sequence of models for H-rich white 
dwarfs (Figures 8 — 10). a) Composition profile showing the mass fractions of C (black), 
O (red), and He (blue), b) Average ion charge (Z), c) Ion coupling parameter T. (Equa- 
tion 8) and electron degeneracy parameter O (Equation 7), d) Ion quantum parameter A 
(Equation 10). 
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models shown, the He layer crosses the warm dense matter regime. Other 
physical parameters, such as the relativity parameter x,, the ion quantum 
parameter A, the total opacity Kio, and the adiabatic temperature gradient 
Vad are very similar to those discussed for the H-atmosphere white dwarf 
models, both in their values and trends, with small differences due to the 
change in mass and charge when He is substituted for H. The ion quantum 
parameter is non-negligible (> 0.01) throughout most of the mass of the 
He layer (only the outer 10~!° of the mass is in the classical limit) and ionic 
quantum effects are apparent in the He convection zone below Tog = 17 000K 
(Figure 12d). The convection zone is associated with the partial ionization 
of He, which corresponds to a maximum of the total opacity and a minimum 
of the adiabatic temperature gradient, just as we observe with the partial 
ionization of H in models of H-rich white dwarfs. 


4.3. Models of hot DQ white dwarfs 


Figure 13 shows the structure of models of hot DQ stars, composed of 
a C-O core overlain with a C envelope (no H or He). Hot DQ stars are 
believed to be the result of the merger of two white dwarfs in a binary system 
and are more massive than typical H-rich and He-rich white dwarfs. The 
models shown are for a mass of 1 Ms. The details of the origin of hot DQs 
and their later evolution are rather uncertain and only models generously 
bracketing their observed temperature range (Zig from 18000 to 24000 K) 
are shown. Little is known about the interior composition of hot DQ white 
dwarfs. In Figures 14 and 15, we adopt a slightly different composition 
where, by construction, the models have equal mass fractions of C and O 
at the center and a decreasing amount of O going outward (Figure 14a). 
Given the low abundance of O seen in the spectra of hot DQ stars, its mass 
fraction is forced to vanish near the surface. Our choice of C-O abundance 
profile is somewhat arbitrary but serves to illustrate significant differences 
with the more common and generally less massive white dwarfs with He 
and H layers. In particular, the higher mass leads to central densities that 
exceed 3 x 107g cm~? where electrons are relativistic. Crystallization has not 
started in these relatively warm models and ions are nearly classical. All 
hot DQ models shown have an outer convection zone associated with the L 
shell ionization and the onset of K shell ionization of C, where the ions are 
strongly coupled but the electrons are non-degenerate. 

The behavior of the average ion charge is more complex than in the H- 
rich and He-rich models (Figure 14b). It starts at (7) = 6.86 at the center 
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Figure 13: Same as Figure 7, but for hot DQ white dwarfs ranging from 15000 to 30000 K. 
The thick black lines indicate where the different electronic shells of C are ionized. The 
envelope is assumed to be made of pure C, and the white dwarf mass is 1 Mo. The surface 
of the star is toward the left, the center is on the right. 


45 


1 Ee as 8 [| teal ies | Tr. t | a a ios | | ri ra 
: : F oma 
20.8 4 el J 
ie) t 4 
Bf 7 [ | 
20.6 jf , +f | 
& ae a + 
L 4 V L 4 
0.4 + L 4 
(00) i =| L 4 
E Ik 1 abt 
al a i 
0 L | @) E | TES BES | | | Te a ER | | eas a | | | ae Te | fl 
0 a =4° -6 - 6. 10 0 5 1G. S15 | =20 
log (1-m(r)/M.) log (1-m(r)/M.) 
Cia at a rAd | Vio Seat free | Ca 8 L! | Cia | Teak | it ae Se | | Ul. eal | 


b 


QO, 


5 [es a) i 
aan yl | = 
E 7 LC 1 
F [ A eer Nae ‘| 
Cuno oe = 
a [ ot See ee ; 
ap OL tra i 
St 4 5 a” 

=6) 4 F 1 
[es } es | | eS = SS: | aS | eS Se: © 4 Lo a ee | as Si | ES SS! | pa} 
0) 5 10 15 20 6) 5 10 15 20 
log (1-m(r)/M.) log (1-m(r)/M.) 


Figure 14: Properties of a cooling sequence of models of a hot DQ white dwarf of M, = 
1 Mo consisting entirely of a C-O plasma. In the absence of a He layer, there is no clearly 
defined core. Models range from Tyg = 30000K to 15000K (heavier line) in steps of 
1000 K. The radial coordinate is a logarithmic mass coordinate (Equation 11) with the 
center of the star located on the left at log g = 0 and the surface towards the right (note 
the “backward” X-axis). In panels b—d, conductive and radiative regions are in blue and 
convectively unstable regions in red. a) Composition profile showing the mass fractions 
of C (black) and O (red), b) Average ion charge, c) Density profile, and d) Temperature 
profile. 
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Figure 15: Same as Figure 14 but showing profiles of the a) Ion coupling (I) and the 
electron degeneracy (Q) parameters, b) Ion quantum parameter A, c) Electron relativity 
parameter x,, and d) Total opacity Ktot. 
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and decreases to (Z) = 6 where the O abundance becomes negligible. The 
ion charge drops to (Z) = 4 as the K shell electrons of C recombine. The 
sharp dip centered at logq = —7 is likely an artifact of the specific C EOS 
table used in these models and hints that a better EOS model for partial 
ionization of C in the WDM regime would be valuable. Finally, most of the 
L shell electrons recombine in the outer region of the models with (Z) = 1 — 
2 at the surface, depending on the value of Teg. The hotter models have two 
distinct convection zones, corresponding to the partial ionization of the K 
and L shell electrons. They are separated by a plateau at (Z) = 4. The two 
convection zones merge when the partial ionization of the two electronic shells 
are no longer distinct in the C EOS, in the (7, p) domain where temperature 
and pressure ionization become comparable. The details of the convection 
zone in hot DQs are closely tied to the modeling of the pressure ionization 
of C in the warm dense matter regime, which can vary considerably between 
EOS tables and has not been benchmarked by experiments. 

The density profiles in the hot DQ models are nearly independent of Tog 
and are about one order of magnitude above those of H-rich and He-rich 
white dwarfs (Figure 14c), a direct consequence of the higher stellar mass. 
The temperature profile is similar to that of H-rich and He-rich models of the 
same Ti range. The merging of the two convection zones at Tog © 18000 K 
is quite apparent in Figure 14d. 

The higher density of the hot DQ models combined with a higher overall 
ion charge (except in the inner core inside logqg = —2) make for an overall 
shift to larger ion coupling (Figure 15a). At the center, I is about 3 times 
larger than in He-rich models for the same Tig and [ > 1 through almost the 
entire star (log g > —11). While the bottom of the convection zone is always 
moderately coupled, it is not degenerate, unlike in the models of H-rich and 
He-rich white dwarfs. The larger density increases the importance of ionic 
quantum effects by about a factor of four (for a given Tyg) and quantum 
effects are non-negligible (A > 0.01) everywhere except for the outer 10~° 
mass fraction of the star (Figure 15b). Similarly, relativistic effects on the 
electrons are now dominant at the center with x, * 2.4 (Figure 15c) and 
remain significant up to the onset of convection. Finally, the profiles of 
the total opacity show the same trends as that of the H-rich and He-rich 
models, except that the two electronic shells of the C atom each produce a 
local maximum in Kio (Figure 15d). The strong correlation between opacity 
maxima and convective instability is also seen in hot DQ models. 

The outer convection zone in white dwarfs plays an important role in the 
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late phases of cooling after convective coupling, in driving the pulsation in- 
stability, and in mixing (Sections 3.3, 3.5). The astrophysics of white dwarfs 
thus depends on a proper modeling of their convective regions, which in turn 
depends, among other factors, on the EOS and opacity of the partially ion- 
ized plasma." Unless the density is very low, partial ionization is challenging 
to model. In H-rich, He-rich and hot DQ white dwarfs, the convective re- 
gion crosses the regime of warm dense matter (Figure 7) whose quantitative 
understanding is actively pursued experimentally and theoretically. 


4.4. Magnetic fields and white dwarf plasmas 

Since about 10% of magnetic white dwarfs have very strong surface fields 
of ~ 10? — 10? MG [22], it is interesting to see what effect, if any, such a field 
has on the interior plasma which is not directly accessible to observation. 

A key quantity in magnetized plasmas is the gyrofrequency of particles 
of mass m and charge Z 

paseo (21) 
mc 
The corresponding energy hw, is the spacing between quantum mechanical 
eigenstates of a free particle in a magnetic field (the Landau levels). Due to 
their much lower mass, electrons have a higher w, and are more affected by 
magnetic fields than ions. 

The effect of magnetization on the EOS of dense plasmas has been studied 
extensively in the context of neutron stars whose typical surface fields are very 
large (~ 10'*G) compared to those found in white dwarfs (e.g. Ref. [117] 
and references therein). The kinetic degrees of freedom of particles in a 
magnetized plasma can be treated classically if ¢ = hw./k,T « 1. Assuming 
B <10°G, this is always true for ions in white dwarfs throughout the entire 
star with ¢ < 10~° in the C/O core and reaching ~ 107° in the very outer part 
of the H layer. The criterion ¢ = 1 is met for electrons only in the very outer 
regions of the star where T’ < 13000K, which falls within the atmosphere 
where the density is very low, the plasma coupling is weak and electrons are 
nearly classical. The corrections to the plasma EOS in the ¢ < 1 regime are 
small but not entirely negligible [167, 117]. Compared to the scale of atomic 
binding energies, a 100 MG field corresponds to an energy of hw, ~ 0.043 Ha. 


While the location of the convection zone boundaries depend on both Vaq and Kp 
through the Schwarzschild criterion (Equations 17-20), the structure and properties of 
the convection zone are largely insensitive to the opacity. 
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Such a field has a small effect on the more tightly bound energy levels of H 
and He, but excited states are moderately affected [168, 169, 170]. The EOS 
of a partially ionized, magnetized H plasma encountered in the atmosphere 
has also been studied [167, 170]. Thus, single white dwarfs with large fields 
fall in a weakly magnetic regime of the EOS where corrections are mostly 
at a perturbative level in the very outer region of the star and are negligible 
throughout the interior. To our knowledge, however, detailed calculations 
of the EOS of H and He in the relatively weak fields encountered in white 
dwarfs have not been applied to models of magnetic white dwarfs. 

On the other hand, highly magnetized white dwarfs with interior fields up 
to 10'*G have been proposed as possible progenitors of rare, over-luminous 
supernovae of Type Ia [171] and as an alternative model to magnetars to 
explain soft gamma-ray repeaters and anomalous X-ray pulsars [172]. A 
simple application of the scalar virial theorem indicates that the global field 
can be as large as 10° G [173] in white dwarfs. Even larger central fields are 
thus physically plausible. Such large fields affect the electron contribution 
to the EOS, which become quantized in terms of the Landau magnetic levels 
[173, 117, 171, 174]. Interestingly, the magnetic field stabilizes the solid phase 
with respect to the liquid and the crystallization temperature increases, with 
the plasma coupling parameter of freezing decreasing from [ = 168 (B = 0) 
to 155 (B = 10'%G) for a carbon plasma [117]. 

With a field of 10'*G, the magnetic pressure contributes significantly 
to the total pressure and affects the structure of the star, notably raising 
the maximum mass of a stable white dwarfs well above the canonical Chan- 
drasekhar mass of 1.40 Mo [175, 164] for a non-magnetic, non-rotating car- 
bon white dwarf [171]. While the Chandrasekhar mass is usually defined in 
terms of the gravitational stability of the star given its EOS (e.g. [176, 175]), 
many factors affect the models near the limit of stability, such as general 
relativity (Section 3.7), the inclusion of the electromagnetic energy in the 
Einstein—Maxwell stress-tensor, the core composition, the equation of state, 
the geometry of the magnetic field, the oblateness of the star and the surface 
boundary condition. At the very high central densities found in stars near 
the maximum mass, instability can arise from two nuclear processes: electron 
capture [176] and pycnonuclear reactions [177]. In models of highly magnetic 
white dwarfs, one of these two nuclear processes usually limits the maximum 
mass. Published estimates of the maximum mass of magnetic white dwarfs 
vary considerably depending on the model assumptions, ranging from 1.75 
to 2.58 Ms, with the more detailed and physically consistent models giving 
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= 2.0 Mo [171, 178, 179, 180, 181, 182, 183, 184, 185, 186], with a correspond- 
ing maximum field of ~ 4 x 10'*G [185]. So far, there is no direct evidence 
for white dwarfs with such high internal fields and there is no known white 
dwarf with a mass above 1.4 Mz [187]. 

Transport properties (diffusion coefficients, conductivity, opacity) are typ- 
ically more sensitive to the details of the microphysics than the EOS. An 
analysis of the physical regimes of transport in magnetized plasmas [188] 
implies that it is unaffected by the magnetic field (unmagnetized regime) 
throughout nearly the entire star. In H-rich white dwarfs, Ht ions become 
weakly magnetized only in the outer ~ 10~'! of the mass of the star for 
fields of 100 MG, and only in stars with Tig 2 20000 K. In the interior, the 
effects of the field are negligible mainly because of the high plasma density 
that results in w, >> we, where w, = \/40Z7e?n/m is the ion plasma fre- 
quency. Thus, ion inter-diffusion remains unaffected (and isotropic) even 
under such strong fields. The field has a stronger influence on electrons, 
which are weakly magnetized over a greater depth in the star, corresponding 
to the outer ~ 10-° M,. Besides the obvious effect on absorption lines at 
the photosphere, the bound-free and free-free opacities are modified at these 
larger depths [189, 190]. We can expect that the screening of the ions by the 
free electrons will become anisotropic and dependent on the field orientation. 
Such a distortion of the electron density surrounding an ion would have an 
indirect effect on the ion inter-diffusion through the effective ion-ion poten- 
tial (Section 5.6.1). While the electron thermal conductivity would also be 
modified slightly by the presence of the field, it no longer plays a role in the 
outer layers of the star. 


4.5. Atmosphere models 


To conclude this section, we examine in more detail the physical regimes 
encountered in the atmospheres of white dwarfs. Despite representing a very 
small fraction of the mass and radius of the white dwarf!?, this thin layer is 
of great importance since it is directly probed by astronomical observations 
via photometry and spectroscopy. Virtually all of the properties of white 
dwarfs are inferred from the analysis of this region. The effects of moderate 
non-ideal physics on the EOS, ionization and dissociation equilibria, atomic 


!2The pressure scale height in white dwarf atmospheres is typically of the order of ten 
to a few hundred meters. 
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level populations and radiative absorption cross-sections are manifest in the 
spectra of white dwarfs, which offer a strong test of microphysics models. 

Figure 16 shows density-temperature profiles for H- and He-atmosphere 
white dwarfs of different surface temperatures. For each profile, a circle in- 
dicates the location of the photosphere. The photosphere is located at a 
Rosseland optical depth of 7, = 2/3, meaning that e~?/? ~ 50% of the pho- 
tons emitted at the photosphere are absorbed or scattered before reaching 
the top of the atmosphere and the other half escapes directly to space. In 
other words, the photosphere is the region from which a typical photon de- 
tected on Earth originates, and it is therefore representative of the conditions 
that we are most interested in modeling correctly. Regions above the pho- 
tosphere also contribute to the emergent flux, particularly at wavelengths 
where the (wavelength-dependent) opacity is high compared to the Rosse- 
land mean opacity (e.g., in spectral lines). The wavelength dependence of 
the opacity implies that the spectrum is formed at different depths in the 
atmosphere. This property can be used to probe the vertical stratification of 
the atmosphere [191]. 

The profiles shown in Figure 16 reveal two important trends: (1) the 
cooler the atmosphere, the denser it is and (2) He atmospheres are denser 
than H atmospheres. To understand why this is the case, we consider the 
definition of the optical depth 7 in terms of the radiative opacity «, dr = 
—kKpdr (the optical depth is measured inward from the surface). The pressure 
in a fluid in hydrostatic equilibrium is determined by 


dP 
= 29 
ip pg (22) 
or 
dP = = dr : (23) 


With the photosphere defined as tT = 2/3, we have, approximately, P = 

29/3. For an ideal gas, the density at the photosphere is of the order of 
2gu 

= ——_.. 24 

 3KkyTe a 

Thus, the photospheric density increases with the surface gravity g and aver- 

age mass per particle ju, and decreases as the photospheric temperature (Ti) 

and opacity increase. The 7g and jz account for part of the trends we see 

in Figure 16 but the opacity plays a major role. Figures 17 and 18 show, for 
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Figure 16: Density-temperature profiles of white dwarf atmospheres. The left set of 
profiles (solid lines) corresponds to H-atmospheres from Teg = 20000 K to 4000 K in steps 
of 1000 K, and the right set (dashed lines) are He-atmospheres with the same Tyg. A surface 
gravity of 10°cms~? is assumed, which corresponds to a typical ~ 0.6 Mo white dwarf. 
The level of the photosphere is indicated by a solid dot. The interior of the star is toward 
the upper right. Note that the profiles shown here differ somewhat from the uppermost 
layers of the profiles displayed in Figures 7 and 11. These profiles were generated using a 
detailed atmosphere code [192, 193] that includes more accurate constitutive physics for 
the specific regime relevant to atmospheres than what is implemented in the more general 
structure code used for Figures 7 and 11. 


53 


two different temperatures, the total wavelength-dependent opacity (black 
lines) of H- and He-atmosphere white dwarfs, respectively. Those figures 
show that (1) for a given atmospheric composition, the atmosphere becomes 
more transparent at low temperatures and (2) for a given temperature, a 
He atmosphere is more transparent than a H atmosphere (note the different 
scales for the opacity in Figures 17 and 18). This is primarily due to the 
higher ionization potential of He, which results in a lower electron density 
and lower opacity. This explains why cool and H-poor atmospheres are more 
dense than warm and H-rich atmospheres. 

Let us now examine the different contributions to the total opacity. In 
the case of a cool H atmosphere (top panel of Figure 17), we see that three 
processes dominate. In the visible, the opacity is dominated by bound-free 
absorption from the H~ ion. In the blue/UV, it is the far red wing of the 
Lyman-a line that contributes the most to the total opacity. Even if it is 
centered in the UV at \ = 0.1215 wm, the Lyman-a line strongly affects the 
emergent spectra of cool H-atmosphere white dwarfs in the visible because 
collisions with H atoms and H2 molecules significantly broaden the absorp- 
tion line profile at those densities. This opacity source has been shown to 
be crucial to reproduce the spectral energy distribution of cool white dwarfs 
[194]. Finally, in the infrared, it is collision-induced absorption (CIA) caused 
by Hy-H» and H»-H interactions that dominates. At sufficiently low tempera- 
tures, CIA becomes so strong that the peak of the spectral energy distribution 
shifts to shorter wavelengths as they cool — the opposite of the behavior of 
a black body spectrum [195, 20]. The physics of CIA is discussed further in 
Section 5.4.2. 

The dominant opacity sources completely change at higher temperature. 
The far red wing of Lyman-a and CIA are both absent from the bottom panel 
of Figure 17. The photospheric density decreases by four orders of magnitude 
between T.¢ = 4000 and 12000 K (Figure 16), which implies fewer collisions 
in the plasma and in turn a less pronounced broadening of the Lyman-a 
line. Similarly, molecular hydrogen is fully dissociated at the lower density 
and higher temperature and Hy CIA goes away. Figure 17 shows that the 
opacity of an H atmosphere at 12000 K is dominated instead by bound-free 
and bound-bound absorption from H. These opacity sources were negligible 
at 4000 K because there was too little thermal energy to excite bound-free 
and bound-bound transitions. Thanks to these new opacity sources, the total 
opacity in the visible is now a few orders of magnitude higher than in the 
4000 K case and the density is lower (Equation 24). Finally, note that the 
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Figure 17: Dominant contributions to the radiative opacity of H-atmosphere white dwarfs. 
The top panel shows the wavelength-dependent opacity of a cool (Teg = 4000 K) white 
dwarf, which would be classified as DC due to the absence of spectral lines, and the bottom 
panel is for a warmer object (Tyg = 12000K), which would be classified as DA due to 
the strong Balmer lines. In both cases, a surface gravity of 10° cms~ 


0.6 


Wavelength (um) 


2 


opacities are evaluated for the conditions at the photosphere. 


ay) 


2.0 


is assumed and the 


He atmosphere, Tog = 4000 K 


| — Total — Het bound-free 
— He scattering — He-He-He CIA 


-- He free-free 


0.6 “10 2.0 
Wavelength (um) 


He atmosphere, Teg = 12000 K 


— Total 
— He bound-bound 


my + : + 


-- He free-free 
— Hey bound-free 
Hey free-free 
10! 
7 
20 
a 
= 
é 
10° 
107! 
0.2 0.3 0.4 0.6 1.0 2.0 


Wavelength (am) 
Figure 18: Similar to Figure 17 but for He atmospheres. The Teg = 4000 K model corre- 


sponds to a DC white dwarf, while the hotter Tg = 12000 K model corresponds to a DB 
(He lines in the visible). 
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opacity of the H spectral lines (grey line) can be a few orders of magnitude 
higher than the continuum opacity set by bound-free transitions (red solid 
line), implying that those lines are formed above the photosphere. 

We now turn to the case of He atmospheres. The most striking aspect 
of Figure 18 is how small the total opacity is for cool He atmospheres. This 
leads to high, liquid-like photospheric densities (Figure 16), where the ideal 
gas physics that is well justified for warmer, less dense atmospheres breaks 
down in the coolest atmospheres with Tyg < 7000K. This complicates the 
modeling of those objects as discussed further in Sections 5.2, 5.3.2, and 5.4. 
For both temperatures shown in Figure 18, the dominant opacity source is 
He~ free-free. This absorption mechanism has a linear dependence on the 
electronic density. This explains why it produces very little opacity at 4000 K 
when He is almost completely neutral, but much more at 12000 K when the 
Saha chemical equilibrium predicts a much higher ionization fraction. Other 
important opacity sources for He atmospheres are Rayleigh scattering from 
He atoms [196] (significant at low temperatures only), CIA from triple-He 
collisions (only at low temperatures, where the density is high enough for 
those three-body collisions to occur), and bound-free and free-free absorption 
from molecular Hej. 


5. Recent advances and challenges in the constitutive physics of 
white dwarfs 


In this section, we describe in greater detail several interesting problems 
where recent progress has been achieved in understanding the physical pro- 
cesses in white dwarfs with computer simulations, better models of the mi- 
crophysics of plasmas, and experimental work that recreates the conditions 
in those stars. The topics cover computational fluid dynamics of convection 
in 3D, equations of state, liquid-solid phase transitions and their phase dia- 
grams, atomic and molecular opacities, line broadening, inter-diffusion, and 
thermal conductivity. 


5.1. Radiation-hydrodynamics simulations of surface convection 


The physics of convective energy transport in stars, including the Sun, is 
relatively well understood [197], although the complexities arising from the 
non-local and turbulent nature of convection has delayed the development of 
precise, quantitative models. Despite presenting a very simple description of 
convection, the mixing-length theory (MLT) has been reasonably successful 
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for white dwarf modelling over the last 50 years. In this picture, a layer is 
either convective or non-convective according to the Schwarzschild criterion 
(Equation 17) or, in the presence of a composition gradient, the Ledoux cri- 
terion [198, 105]. In reality, the convective flows do not stop abruptly at 
well-defined boundaries at the top and bottom of the convection zone. Fur- 
thermore, in the non-adiabatic convective layers that define the atmosphere 
of most stars and white dwarfs, the convective efficiency predicted by the 
MLT is very sensitive to the underlying parameters describing the geometri- 
cal shape, radiative energy losses, and lifetime of convective cells [199, 130]. 
These quantities are not constrained by the MLT and must be calibrated 
from observations. 

The emergence of powerful computing resources has enabled the gen- 
eration of large grids of 3D radiation-hydrodynamics (RHD) simulations 
that have changed the paradigm for convection in the outer layers of stars 
[200, 201, 202]. Large grids of 3D simulations of convection are now covering 
almost all branches of stellar evolution including main sequence stars, brown 
dwarfs, giants, supergiants and white dwarfs [203, 204, 135, 205], enabling 
studies of large samples of stellar spectra with 3D atmosphere models. Two 
independent state-of-the-art codes have produced the bulk of the 3D RHD 
models for white dwarfs: CO5BOLD [203] and ANTARES [206]. 

These RHD simulations solve the time-dependent fluid equations of con- 
servation of mass, momentum, and energy in three spatial dimensions for 
a compressible plasma in a constant gravitational field. This is solved to- 
gether with non-local, frequency-dependent radiation transport and tensor 
viscosity using a realistic EOS. The radiative opacities can be either gray (i.e. 
frequency averaged) or chosen to preserve some level of frequency-dependent 
via an opacity-binning scheme. For both main sequence stars and white 
dwarfs, the box-in-a-star simulation geometry is used where the box repre- 
sents a small portion of the surface of the star. A Cartesian grid is setup 
with periodic lateral boundary conditions and appropriate physical boundary 
conditions at the top and bottom (e.g. constant radiative or entropy flux). 
These simulations have been compared to observations of the surface of the 
Sun where resolved convective eddies are apparent as granulated brightness 
fluctuations at optical wavelength. They have been shown to agree within a 
few percent in brightness contrast [207]. 

Surface convection in white dwarfs (Figure 19) is quantitatively similar to 
that seen in main sequence stars including the Sun since the Mach (0.01-1.0) 
and Péclet (10~*-10*) numbers as well as the granulation intensity contrasts 
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(0.1-25%) are much alike [208], although with some notable differences for 
He atmospheres [205]. The major difference between white dwarfs and main 
sequence stars is in the extent of the internal convective envelope. The radial 
extent and total mass enclosed in white dwarf convection zones are orders of 
magnitude smaller [130]. Furthermore, the higher pressure at the surface of a 
white dwarf, combined with the usually pure-H or mixed H/He atmospheric 
compositions, results in different 3D effects on the spectral line shapes [135]. 

Overall, these 3D simulations have demonstrated that it is possible to 
determine largely parameter-free, time- and spatially-averaged temperature 
and pressure stratifications for convection zones, in contrast to the 1D MLT 
theory where free parameters must be fitted to the observations. Typically, 
no single 1D MLT parametrization can reproduce the full surface-averaged 
emergent spectrum of a 3D RHD model of a white dwarf at all wavelengths 
(Figure 20). The reasons for this behaviour are complex and depend on how 
3D fluctuations (e.g. cool downflows and warm updrafts) average over time, 
the stellar disk, and optical depth. 

Current grids of 3D model atmospheres cover essentially the full range 
of convective white dwarfs, from 3500 to 18000 K and log g from 5.0 to 9.0 
for pure-H atmospheres, and from 10000 K to 34000K and logg from 7.5 
to 9.0 for mixed H/He atmospheres [210, 205]. Only the warmest pure- 
H (211000K) and mixed atmospheres (224000 K) currently have the full 
vertical extent of the convection zone simulated, allowing the study of phys- 
ical processes at the bottom boundary layers, such as convective overshoot 
[206, 129] and pulsation excitation mechanisms. With current computing 
capabilities, simulations of the full extent of the convection zone of cooler 
white dwarfs remain out of reach. 

Accounting for 3D convection in calculations of stellar structures is de- 
sirable for many applications, such as stellar evolution and asteroseismology. 
The calibration of 1D envelope structures with 3D simulations as surface 
boundary condition for interior models has now been performed for white 
dwarfs [130, 205] using methods previously employed for main sequence stars 
[199, 211]. For warmer convective white dwarfs, the full vertical extent of 
the convection zone can be simulated in 3D, which can then be matched 
directly with the appropriate 1D structure below the convection zone. For 
cooler white dwarfs, the 3D simulations of the upper part of the convection 
zone can be matched to 1D MLT models of convection to cover its full ex- 
tent. Deep convection becomes essentially adiabatic, i.e. constant specific 
entropy with depth, a limit whose thermodynamic stratification can be well- 
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Figure 19: Snapshot of a 3D radiation-hydrodynamics simulation of convection in a white 
dwarf at Tog ~ 12000K and logg = 8. Top: Temperature structure for a slice in the 
horizontal-vertical 2 — z plane through a rectangular box with coordinates x, y, z (in km), 
with z being in the vertical direction. The temperature is colour coded from 60 000 (red) to 
7000 K (blue). The arrows represent convective velocities, while white lines correspond to 
contours of constant Rosseland optical depth 7,, with values given in the figure. Bottom: 
Frequency integrated emergent intensity at the top of the horizontal « — y plane. This 
pattern of bright cells with dark borders, called granulation, is observed at the photosphere 
of the Sun. Hot rising convective flows form the bright areas, where the fluid cools by 
radiation and sinks along the darker edges. The rms intensity contrast with respect to the 
mean intensity is 18.8%. The length of the bar in the top right is 10 times the pressure 
scale height at the photosphere (7p = 2/3). Reproduced from Ref. [209] with permission 
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Figure 20: Comparison of the blue wing of six Balmer line profiles (Ha to H8) calculated 
from spatially and time-averaged 3D atmosphere structures (red, solid) and standard 1D 
structures (black, dashed). The values of Teg are shown on the different panels and 
log g = 8. All line profiles were normalized and shifted vertically for clarity. The overall 
trend is that the 3D RHD simulations produce lower Balmer lines that are narrower, and 
the opposite for the higher lines in the series. Reproduced from Ref. [209] with permission 
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represented with the MLT. The constant value of the entropy in the regime 
of adiabatic convection can be extracted at the bottom of a 3D RHD simu- 
lation to calibrate the 1D MLT that is used to extend the convection zone 
along an adiabat all the way to the bottom. Figure 21 demonstrates that 
significant entropy fluctuations are apparent at all geometrical depths in a 
3D simulation, although there is a constant asymptotic upper limit below 
the observable surface. According to Ref. [212], the fluid in central regions 
of broad ascending flows is still thermally isolated from its surroundings as it 
rises adiabatically until it reaches layers immediately below the photosphere. 
In other words, convective upflows keep an imprint of the physical conditions 
at the bottom of the convection zone. The asymptotic entropy value of the 
3D simulations has therefore been employed to calibrate 1D structures in 
established white dwarf evolution codes. 


5.1.1. Convective overshoot 

The size of a 1D convection zone is based on the location of the Schwarzschild 
instability boundaries. It is well-known that this approach cannot properly 
characterize the actual hydrodynamic flows that exist at convective bound- 
aries, instead producing a significant and unphysical discontinuity in the flow 
of kinetic energy [213, 214, 201]. 

In reality, plumes are accelerated near the edge of the unstable layers 
and are able to travel outside of the 1D convection zone into the convec- 
tively stable layers in a phenomenon known as convective overshoot, which 
happens both at the top of the convection zone, usually in the photosphere, 
and at its base. A robust characterization of convective overshoot requires 
hydrodynamics simulations. The strength or characteristic length of convec- 
tive overshoot is usually different when looking at thermal (convective flux) 
or mechanical quantities (convective velocities). One useful definition of the 
extent of overshoot corresponds to the region where the macroscopic over- 
shoot diffusivity is larger than the microscopic diffusion coefficient. These 
regions are considered to be fully mixed in terms of chemical abundances, 
from which follows the definition of the mixed convective mass. Recent inde- 
pendent studies [206, 129] of warm H-rich white dwarfs with small convective 
zones have demonstrated that bottom convective overshoot has a dramatic 
impact as it increases the inferred mixed mass by up to 2-3 orders of mag- 
nitude (see Figure 22). Bottom convective overshoot has strong implications 
for spectral evolution (convective mixing) and accretion of planetary debris 
(inferred accretion rates and mass of accreted planetesimals). The impact 
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Figure 21: Local 3D entropy values (black dots) as a function of geometrical depth for 
a subset of a simulation at Teg = 10025K and logg = 8. The entropy profile of the 3D 
simulation, averaged over constant geometrical depth z and time, is shown with a red 
solid line. The entropy of an equivalent 1D model with the best MLT parameterization 
is shown as a red dashed line. Average optical depth values of Tt; = 100, 1, and 0.1 are 
highlighted with dashed blue lines. The surface of the star is to the right. Adapted from 
Ref. [130], with permission @) AAS. 
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of convective overshoot on the convective coupling to the degenerate core 
at lower Tyg has not been explored but is likely small compared to existing 
uncertainties in the constitutive atomic and equation of state physics. 

Convective overshoot studies have so far been limited to white dwarfs with 
Teg = 11000 — 18000 K and shallow convection zones. There is nevertheless 
evidence that overshoot becomes less important in terms of the increase in 
mixed mass as the convection zone reaches deeper layers [129] but more work 
is needed. A physically motivated analytical parameterization of convective 
overshoot [144] is of major importance for white dwarfs and stellar evolution 
in general. 

In cases where material is accreted onto a white dwarf from rocky plan- 
etesimals or a stellar companion, a gradient in chemical composition is in- 
duced between the photosphere (or mixed layers) and the envelope just be- 
low. This gradient can lead to thermohaline instabilities [215, 216, 144]. This 
mechanism is distinct from convective overshoot, and can happen in warm 
white dwarfs that have no convective instabilities. While no 3D simulations 
have yet been performed for thermohaline mixing in white dwarfs, recent 
studies have shown that it is likely to be a significant effect adding up to 
convective overshoot [144]. 


5.1.2. Magnetic white dwarfs 

It has long been predicted—and observed—that magnetic regions in the 
Sun and low-mass stars locally inhibit surface convection and convective en- 
ergy transport [217, 218, 219]. Based on a handful of magneto-hydrodynamics 
(MHD) calculations for white dwarfs at Tg ~ 10000 K and pure-H compo- 
sition, it has been shown that magnetic fields as low as 1kG can inhibit 
surface convective energy transfer [137]. This is illustrated by the disappear- 
ance of the fluctuations in surface intensity (granulation) as the field strength 
increases (Figure 23). This results in a steeper vertical temperature gradi- 
ent in the atmosphere, which can be constrained using spectroscopy [220]. 
The external magnetic fields necessary to inhibit convection at the surface 
varies from 1kG for warmer DA white dwarfs to almost 1 MG for cool He- 
dominated atmospheres [221]. It is important to emphasize that inhibition of 
bulk convective energy transport does not mean inhibition of small scale con- 
vective motions and instabilities. Three-dimensional MHD simulations show 
that inefficient convective flows act on smaller spatial scales in the presence 
of a magnetic field, but still show large local convective velocities [137]. This 
suggests that convective chemical mixing is still effective even for magnetic 
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Figure 22: Diffusivity as a function of depth log gq (Equation 11) in a pure-H white dwarf 
envelope with Tag = 12009K and logg = 8. The macroscopic (convective) diffusivity 
inferred from 3D convective velocities [129] is shown with solid lines. From top to bottom, 
the diffusivity is shown with increasing spatial filtering of the velocities to remove small 
scale features. This is done to reduce the effect of artificial standing waves. The best fit 
to macroscopic diffusion coefficients is shown as a black dot-dashed line. The convectively 
unstable region defined by the Schwarzschild criterion is indicated by the orange shaded 
region, which has a much smaller extent than the high-diffusivity profile extracted from the 
simulation that includes convective overshoot. Also shown are the microscopic diffusion 
coefficients for Mg, Ca, and Fe (dotted lines), taken as typical trace elements in white 
dwarf atmospheres. The base of the mixed region is located where the macroscopic and 
microscopic diffusion coefficients are equal (filled color points in 3D and black point in 
1D). The surface of the star is to the left. Reproduced from Ref. [129] with permission. 
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fields that are one or two orders of magnitude stronger than necessary to 
inhibit convective energy transfer. This has implications for the visibility of 
metals in the study of metal-polluted atmospheres. 


5.1.8. Horizontal mixing and rotation 
In addition to efficient vertical mixing of chemical elements, convection 
can also mix material in the horizontal (or surface plane) direction. This has 
received little attention until recently [222, 221], in part because this process 
is bound to be relatively fast given surface convective velocities reaching 
~ 1kms~! (both in the horizontal and vertical directions), resulting in a 
horizontal diffusion time scale of the order of 
R? 7 R? 


tspread ™ 
P D 


conv Disa Ls ~ : - , ey 
where Dyeony is the convective diffusivity and H, the pressure scale height. 
The most important parameter is the ratio between the horizontal and verti- 
cal diffusion time scales, which describes whether accreted material has time 
to spread horizontally before it sinks and disappears from the visible layers. 
It has been found that most white dwarfs are able to efficiently spread ac- 
creted material, with the exception of DA white dwarfs above 11000 K. The 
addition of thermohaline mixing could make horizontal spreading somewhat 
more efficient, making that transition a few thousand degrees higher. 

This prediction is possibly in tension with observations of warm (= 11000 K) 
polluted DAZ white dwarfs, which show strong evidence of a homogeneous 
metal-polluted surface [221]. The solutions to this problem remain open, 
but depend on the unknown geometry of accretion. In particular, near- 
spherical deposition would remove the need for an unknown source of hori- 
zontal spreading to explain the observations. Magnetic white dwarfs accret- 
ing planetary debris will be of particular interest in that context, allowing one 
to constrain both accretion geometry and mixing in the presence of magnetic 
fields. 

Rotation is slow in single white dwarfs. Ref. [223] finds an average pe- 
riod of 35h with a standard deviation of 28h from studies of pulsating white 
dwarfs. Similar results are found using periodic variability of spotted mag- 
netic white dwarfs [224]. Such slow rotation is not known to have observable 
effects on the structure, cooling age, surface composition, or magnetic field 
generation (except for the recent scenario of a crystallization dynamo [159]). 
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Figure 23: Frequency-—integrated radiative intensity emerging from the x — y horizontal 
plane at the top of the computational domain for 3D MHD simulations. All simulations 
have a constant surface gravity of log g = 8.0, a pure-H composition, and rely on the same 
entropy value for the inflowing material through the open bottom boundary, resulting in 
Teg ~ 10000 K. The panels show simulations where a constant vertical field of 0, 0.5 and 
5kG was imposed at the lower boundary. The rms intensity contrast with respect to the 
mean intensities is also shown above the panels. The length of the bar in the top right is 10 
times the pressure scale height at 7, = 2/3. Reproduced from Ref. [137] with permission 
@©AAS. 


67 


5.2. Non-ideal chemistry in cool He-dominated atmospheres 


The liquid-like density (0.1 — 1gcm~?) that characterizes cool He-rich 
atmospheres (Figure 16) poses a significant challenge to their modeling. The 
equation of state P(p, 7) of He is already satisfactorily described by existing 
tables [225]. The key problem is instead characterizing the chemical equi- 
librium of minor ionic and molecular species that do not affect the equation 
of state but contribute disproportionately to the opacity under those con- 
ditions. The challenge is to calculate, for a dense interacting He fluid at 
density p and temperature T, the ionization ratios (e.g., n(He*)/n(He)) and 
the molecular dissociation ratios (e.g., n(H)/n(Hz2)) of the different atomic 
and molecular species present in the atmosphere. This is a problem of the 
utmost importance for the study of those objects, which include some of 
the most intriguing white dwarfs, such as the IR-faint (or ultracool) white 
dwarfs and the recently discovered Li-polluted white dwarfs. The chemical 
equilibrium prescription used in an atmosphere model calculation determines 
the radiative opacities (e.g., more Hz results in more CIA), and in turn the 
structure of the atmosphere and the resulting spectrum. 


5.2.1. The chemical picture 

All stellar atmosphere codes currently rely on the chemical picture. In this 
approximation, bound species (atoms and molecules) are considered as basic 
particles. This is to be contrasted with the physical picture, where only nuclei 
and electrons are the basic particles. Assuming that atoms and molecules are 
well-defined entities can be problematic at high densities where the boundary 
between bound and free states is blurred. The abundance of a molecule 
can depend on the very definition of a “molecule” [226]. In spite of this 
fundamental limitation, the chemical picture is much more practical when 
it comes to model atmosphere calculations. To include atomic spectral lines 
and molecular bands, the population of thousands of atomic and molecular 
levels is needed, making the physical picture impracticable. 

Under low-density conditions, the ionization and dissociation ratios are 
obtained with the Saha equation, which is derived assuming an ideal gas 
chemical potential for all species involved, 


Ni; h? 


Q; (QamikpT)?” 


we = Exo + kT In ; (26) 


with Ej; the ground-state energy of species i and Q; its internal partition 
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function, 
Eig 
Q;= i Gi,j EXP (-=) ; (27) 


where the sum runs over all bound states of species i, and where g;,; is the 
statistical weight of state j.1° 

The chemical equilibrium is obtained by equating the chemical potential 
of reacting species. For example, for the reaction H = Ht + e, we have 
Ly = [yt + Ue. The system of equation is closed by imposing overall charge 
neutrality and the conservation of the number of atoms of each element. At 
the densities encountered in cool He-atmosphere white dwarfs, two modifi- 
cations must be introduced: (1) an interaction (or non-ideal) term has to 
be added to Equation 26, and (2) the partition function should be modified 
to account for the perturbation of the internal degrees of freedom due to 
interactions with other particles. As much as possible, the latter should be 
consistent with the former. 


5.2.2. The occupation probability formalism 

During the last three decades, the occupation probability formalism of 
Hummer & Mihalas [227] has been widely used to evaluate chemical equilib- 
rium under non-ideal conditions in white dwarf atmospheres [228, 134, 229, 
230, 231, 100]. In this framework, an occupation probability w; is assigned 
to each bound state with w; = 1 for unperturbed levels and w; = 0 for levels 
that are completely destroyed by interactions. Those w; then multiply each 
term in the Boltzmann distribution, 


a B 


NY 


and in the calculation of the partition function, 


Qi = 2 Wisdhs exp (-=4) . (29) 


13 lectrons are normally classical in the atmospheres of white dwarfs but Equation 26 
can be easily modified to account for degeneracy (see Figures 7, 11, and 13). 
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For cool He atmospheres, the w;,; are dominated by neutral interactions with 
neighboring particles. In this case, Hummer & Mihalas give the prescription 


Wij = EXP -= Song (rng Hee) | (30) 
vig! 

where r;,; is a radius assigned to species 7 in state 7, and the sum runs over 
all species and all states.‘ While very simple to implement and interpret, 
the excluded volume effect of Equation 30 is only a crude approximation of 
the actual interaction between two neutral particles. Moreover, there is no 
theoretical prescription for the radii r;,;; they are free parameters. The r;,; 
can be arbitrarily adjusted to obtain a good agreement between the white 
dwarf spectra and model atmospheres [228]. However, because of the lack of 
predictive character of this approach, a good agreement between theory and 
observation does not mean that the atmosphere is being properly modeled 
nor that the derived atmospheric parameters are accurate. 


5.2.3. More advanced interaction models 
Recently, chemical equilibrium calculations that are more predictive than 
the Hummer & Mihalas formalism have been performed [233, 234, 235, 192]. 
Those works have focused on accurately computing the non-ideal chemical 
potential p24 that has to be added to Equation 26 to account for interactions 
in a dense He fluid, 
nid 


bse = we + ue (31) 
In general, p24 


nid is a function of JT and n;. Adding those p?'¢ to each species 
i involved in the ionization (dissociation) process being studied effectively 
shifts the ionization (dissociation) potential. There are two contributions to 
consider in the calculation of uw". The first is simply the interaction energy 
of species 7 with the He fluid, and the second is an entropy term. The latter 
arises because the surrounding medium responds to the ionization or dissoci- 
ation of a particle: the radial distribution function g(r) of the fluid depends 
on the ionization or dissociation state, which in turn affects the chemical 
potential. Both contributions can be evaluated using the classical theory of 
fluids [236, 237], provided that realistic pair interaction potentials are avail- 
able. Such calculations have been performed for the chemical equilibrium 


l4The argument of the exponential corresponds to the linear contribution in the packing 
fraction of the hard sphere equation of state [232]. 
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of H/H, [233], H/H~ [234], Het/He [235], and Ct/C, Nat/Na, Mgt/Mg, 
Cat/Ca, and Fet/Fe [192] in a bath of He. 

Despite this recent progress, many aspects of current chemical equilib- 
rium models should be improved. A first limitation of current calculations 
is that some are based on ab initio pair potentials obtained in the infinite 
dilution limit. This is problematic as this approach does not account for the 
many-body effects that become important at the liquid-like densities con- 
sidered here. This problem can be solved by using potentials calibrated to 
high-pressure experimental data [235] or by using Density Functional The- 
ory Molecular Dynamics (DFT-MD) to simulate the behavior of a many-body 
system at the relevant densities [192]. 

A second limitation with the current state of the art is that the inter- 
nal partition functions Q; are assumed to remain unaffected by the inter- 
actions. All calculations are performed assuming that all species remain in 
their ground states, and, unlike the occupation probability formalism, noth- 
ing is said about differences in how the different internal energy levels of 
each species are affected by non-ideal effects. This also means that there is 
currently no accurate prescription for how the Boltzmann distribution (Equa- 
tion 28) is perturbed by non-ideal effects. Knowing the population of each 
internal energy level of a species is needed to predict the strength of spectral 
absorption lines and bands, making this problem particularly important. 


5.2.4. Pressure ionization and the phase diagram of helium 

Pure He atmosphere models have an interesting history related to the 
difficulty of properly modeling the pressure ionization of He at low tempera- 
tures (T’ < 1eV), where the ionization fraction of He is very small (< 107°) 
but highly consequential and very sensitive to the details of the ionization 
model. Atomic He has a very low opacity at optical and infrared wavelengths 
(a combination of He Rayleigh scattering and He-He-He collision-induced 
absorption, Figure 18) resulting in high-density atmospheres. The presence 
of a small amount of free electrons, from the ionization of He or from traces 
of other elements with relatively low ionization potentials, introduces He~ 
free-free electron scattering and He} bound-free absorption that increase the 
opacity considerably and decrease the overall density of the atmosphere (Fig- 
ure 18). 

Early models [238] used an EOS and opacities based on the finite-temperature 
Thomas—Fermi—Dirac model of the isolated atom for the intra-atomic poten- 
tial. In retrospect, this is a rather crude model for this (p, 7) regime, lead- 
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ing to a poor EOS and exotic atmosphere models where conductive energy 
transport is significant at the photosphere. A new set of pure He atmosphere 
models [134] introduced He-He interactions in the calculation of the chemi- 
cal equilibrium using a He-He pair potential calibrated to shock data up to 
56 GPa [239]. The latter ensured a realistic P(p,T) relation for dense He 
that also agreed quite well with ab initio simulations of dense He [235] up 
to p+ 0.8gcem-3. This allows the evaluation of 4 (Equation 31), which 
shifts the ionization equilibrium with respect to the ideal Saha equilibrium. 
The He-He interaction being strongly repulsive, wt > 0 and the ionization 
fraction of He increases. The new models also included He Rayleigh scatter- 
ing and He~ free-free opacities. The resulting low-7.,g¢ models had densities 
as high as 4gcm~? and pressures of 20TPa at temperatures below leV. 

Based on those high densities and pressures, the intriguing possibility 
that phase transitions in He at high pressures occur in white dwarfs was ex- 
plored. Two recent ab initio calculations pertain to this question. One is the 
insulator-metal transition of solid He [240], which occurs around 35 TPa, the 
other is the melting curve of atomic He [241]. However, the improved ioniza- 
tion model of Ref. [235], which also includes He-e and He-He* interactions 
— thus providing a more complete description of the ionization equilibrium 
of He by introducing additional contributions to u%? and additional charged 
species — predicts an increase of the ionization fraction and of the He~ free- 
free opacity of about 3 orders of magnitude. Consequently, current pure He 
atmosphere models have maximum densities of ~ 1 gcm~°? and barely reach 
100 GPa. As a result, the models do not approach the melting and metal- 
lization coexistence curves and remain well into the fluid/plasma regime; no 
phase transition is expected (Figure 24). 

The sensitivity of the atmospheric profile to the very small ionization 
fraction illustrates the importance of accurately modeling the onset of pres- 
sure ionization of He in cool atmospheres. We stress that the ionization 
model of Ref. [235], which represents the state of the art in the context of 
white dwarf atmosphere models, remains provisional. An indirect validation 
of the ionization model of He might be possible by measuring the opacity 
of He compressed dynamically or in a laser-heated diamond anvil cell (Sec- 
tion 5.8.3). 

A third set of simulations focused on the transition of He from an insu- 
lating atomic fluid to a conducting, fully ionized plasma [242]. This recent 
work found that the EOS, the conductivity and the ionization fraction all 
vary smoothly as density is increased at fixed temperature. A useful defini- 
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Figure 24: Comparison of recent calculations of the phase diagram of He with white dwarf 
models. The high pressure melting curve of He [241] (heavy blue curve) and the insulator- 
metal transition in the solid phase [240] (heavy purple curve) are first order transitions. 
The insulator-metal transition in the fluid phase (also known as pressure ionization) is 
indicated by the closure of the band gap (blue dotted curve) and occurs continuously 
[242]. These transitions in dense He are compared to pure He white dwarf atmosphere 
models of Tag = 4000K to 8000K (left panel, orange lines) and He-rich white dwarf 
interior models with Teg = 4000 K and 5000 K (right panel, red curves). All models shown 
have a typical white dwarf surface gravity of 10° cms~?. The models stay clear of the first 
order phase transitions but those with Teg < 6000 K cross the regime of pressure ionization 
(right panel) at depths well below the atmosphere. 1TPa = 10 Mbar = 10% ergcm™?. 


tion of pressure ionization that is borrowed from condensed matter physics 
is the density at which the band gap between the valence and conduction 
electron bands closes, which corresponds to a vanishing ionization energy. 
Figure 24 shows the calculated locus of (p, 7’) points of the band gap closure 
of He. As we saw in Figure 11, the coolest He-rich white dwarf models cross 
the regime of pressure ionization but at a depth well below that of the atmo- 
sphere. This careful study of pressure ionization provides a benchmark for 
tables of the EOS and transport properties of He in the regime of 1-20 gcm™3 
and 1—5eV. Note that the method used is not suited for quantifying the very 
low ionization fraction in cool He-rich atmospheres. 
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5.3. Theoretical atomic line profiles 

The determination of the basic physical characteristics of white dwarfs, 
such as their surface temperature, gravity, and composition, from which their 
mass, radius and age are derived, is usually based on a detailed analysis of 
their spectra. The method of fitting the absorption lines developed three 
decades ago [243] has become a work horse in the field and has been applied 
to over 18000 white dwarfs (Figure 25). The success of this spectroscopic 
method depends on the reliability of the modeled line profiles, primarily in 
the UV and optical. This puts stringent demands on the modeling of pressure 
broadening, particularly of H and He lines, the two most common elements 
seen in white dwarf spectra, and on the atomic physics of dilute and weakly 
non-ideal plasmas. In the fully or partially ionized plasmas typical of white 
dwarf atmospheres, Stark broadening dominates but in the coolest white 
dwarfs broadening by collisions with atoms (H or He) and even molecules (H2) 
takes precedence. Because of the relatively simple atomic structure of H and 
He, they have been studied extensively both theoretically and experimentally. 

Note that natural and thermal (Doppler) broadening are generally negligi- 
ble in white dwarf atmospheres and pressure broadening dominates. Natural 
line broadening is due to the uncertainty principle, AE ~ h, where T is the 
atomic state lifetime. Typically, we have tT ~ 1078s and 

AE Av h 


= oon ae y) 
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for a transition at 5000 A. Thermal broadening is given by 


A erma 2kpl x 
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for aH gas at 20000 K. From those estimates, we see that natural and thermal 
broadening do not broaden spectral lines by more than 1 A, a much smaller 
value than the line widths observed in white dwarfs (Figure 25). 


5.8.1. Broadening by collisions with charged particles 

Interactions between a radiating atom or ion with ions and electrons in 
the surrounding plasma lead to Stark broadening of spectral lines. Even 
after decades of studies of atomic line broadening in plasmas, the theory, 
as implemented in several line shape codes, remains unsatisfactory and a 
work in progress [244]. The theory is complex, which motivates several ap- 
proximations of unknown accuracy that are variously implemented in codes 
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Figure 25: Example of the spectroscopic method applied to fitting the Balmer lines of 
the DA white dwarf WD 0307+149. Data are shown in black and the best fitting model 
(with Tog = 22285K and log g(cms~?) = 8.04) in red. The five individual absorption 
lines shown have been extracted from a continuous optical spectrum. Figure courtesy of 
P. Bergeron. 
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[245, 246]. This is of concern as Stark broadening is used as a diagnostic of 
conditions in laboratory plasmas and various codes/models lead to different 
inferred values of the electron density. The problem is compounded by the 
very limited amount of experimental data that can benchmark the theory 
and even fewer that have independent diagnostics of the temperature and 
electron density (Section 5.8.1). 

Applications of the unified Stark broadening theory [247] gives good fits 
of the spectra of DA white dwarfs (Figure 25) but fits of individual Balmer 
lines leads to systematic differences in the derived stellar parameters. Closer 
attention to the perturbed population of the H atom levels in a slightly non- 
ideal gas regime partly resolved this internal inconsistency [231], yet resulted 
in a white dwarf mass distribution that is dependent on Teg, which is unrea- 
sonable on astrophysical grounds [51]. Modeling the line formation through 
an inhomogeneous convective atmosphere (Section 5.1), rather than under 
the standard assumption of a one-dimensional, vertically stratified atmo- 
sphere model, appears to have largely resolved this issue [135]. Nonetheless, 
this sparked a renewed interest in the experimental measurement and the 
theory of line profiles of H and He relevant to white dwarfs that revealed 
some surprises. 

Recent developments in the theory of Stark broadening with a view to- 
ward applications to denser plasmas and white dwarfs started with an investi- 
gation of the common dipole interaction approximation [248]. The Coulomb 
potential energy between a perturber particle (ion or electron) and the ra- 
diator’s nucleus and electrons can be expanded as a series of multipole in- 
teraction terms, which is usually truncated after the second (dipole) term. 
By evaluating higher order terms and treating the ions and electrons on the 
same footing, it was found that the dipole approximation is inaccurate and 
that the quadrupole term is important to match the experimental line profile 
of HB at n. ~ 10'8cm~? and T = 1 eV [248]. The line profile is converged to 
< 1% when including the octupole term. When using H( as a diagnostic of 
the electron density under those conditions, the quadrupole moment brings 
a 12% correction to ne. 

A series of improvements followed. The application of a different numer- 
ical method to solve the Hartree-Fock equations improved the accuracy of 
the wave functions of the bound and the perturber electrons used to calcu- 
late the transition probability [249]. The inclusion of correlations between 
the radiator and the plasma in the relaxation theory of electron broadening 
[250] resulted in a considerable improvement in the agreement with the line 
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shift for ne = 2-10 x 10'°cm~%. This effort has culminated in a systematic 
study [246] of the three main approximations in Stark line shape model- 
ing: 1) The second order expansion of the broadening operator (T-matrix), 
2) The dipole approximation of the radiator-plasma interaction (see above), 
and 3) The classical treatment of the perturber electron. Various combina- 
tions of these approximations have been used in previous work but this is 
the first to have lifted all three. Besides allowing a systematic assessment 
of the accuracy of each approximation, the application of this new model 
to the Lyman a and Lyman £ lines in a hydrogen plasma (T = 1-7eV, 
ne = 10'8cm7?) shows significant differences with the model widely used in 
white dwarf atmospheres [247] and the more recent work removing the second 
order approximation [248]. These recent improvements in Stark line shape 
theory have proved very valuable in the analysis of the He@ line in a highly 
ionized, solid density Ti plasma (Ton © 600eV, n. © 1.1 x 10%4cm~?) [251]. 
This systematic progress in the theory and calculations of Stark line profiles 
in weakly coupled plasmas is an important development in plasma physics". 
A preliminary implementation of the latest line profiles for Ha through Hy 
shows very little change in the synthetic spectra of DA white dwarfs but they 
are found to be sensitive to the treatment of the occupation probability and 
its effect on the pseudo-continuum opacity [253]. The impact on the analysis 
of white dwarf spectra has not yet been assessed. 

The more complex case of He lines in the optical was recently revisited 
[254]. For over two decades, the Stark-broadened He line profiles used in 
white dwarf atmosphere codes have been those of Ref. [255], who applied 
the semi-classical (and semi-analytical) theory of Stark broadening to 21 lines 
in the optical. This standard approach to Stark broadening relies, among 
other approximations, on the quasi-static approximation for the ions and 
the impact approximation for the electrons. The latter was complemented 
with the so-called ”one electron” approximation to improve the transition 
between the impact approximation for the electrons, which is valid in the 
line wings, and the quasi-static approximation for the electrons in the line 
core. By applying a simulation method to model the collisions between the 
radiating atom and the perturbers (ions and electrons, treated as classical 
particles), it becomes possible to do away with the quasi-static approximation 
for the ions and include the dynamics of ion-radiator collisions on the same 


15Ref. [252] gives a comprehensive introduction to the current state of the theory. 
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footing as electron-radiator collisions. Weak and strong collisions are also 
treated without distinction, resulting in a naturally self-consistent description 
of the line core and line wings [254]. This simulation approach is similar 
to that used for modeling Stark broadening of H lines [253, 248] although 
Ref. [254] retained many common approximations in the theory of Stark line 
broadening, such as the dipole interaction between radiator and perturber. 
Six optical lines of neutral He were modeled under a range of conditions found 
in white dwarf atmospheres. As was found for the lines of H, the ion dynamics 
lead to variations from the profiles of Ref. [255], particularly near the line 
core and the forbidden components, but the effect is very small. Inclusion 
of the new profiles for the 4471 A and 4922 A He lines in atmosphere models 
shows minimal deviations from models spectra based on the standard profiles. 
It will be interesting to see whether increasing the level of physics fidelity 
by lifting other important approximations — as was done for hydrogen — 
will lead to significant effects in He lines. As is the case for the new H Stark 
line profiles, additional calculations are needed to quantify the impact of this 
work on the analysis of white dwarf spectra. 


5.8.2. Broadening by collisions with neutral particles: Cool white dwarfs 

In cool stellar atmospheres, the broadening of atomic absorption lines is 
dominated by collisions between the radiating atom and perturbing atoms 
or molecules (predominantly H, He, and Hz). This broadening mechanism 
can usually be described using a simple Lorentzian profile. This assumption 
is well justified in the low-density limit (9 < 10~* gcm~?), where the impact 
approximation can be made, i.e., when the duration of a collisions is small 
compared to the time elapsed between collisions. However, cool white dwarf 
atmospheres (especially those dominated by He) can reach much higher den- 
sities (Figure 16). Under those conditions, the impact approximation breaks 
down, and a more sophisticated approach is required to describe the colli- 
sionally broadened lines. In this regime, white dwarfs exhibit severely dis- 
torted and shifted absorption lines, and Lorentzian profiles are a very poor 
description of the observations. This is particularly problematic for cool DZ 
white dwarfs, where an accurate modeling of the spectral lines of the heav- 
ier elements is needed to reliably extract the composition of the accreted 
exoplanetary material. 

The so-called unified line broadening theory has been developed to de- 
scribe collisionally broadened line profiles under those challenging conditions. 
It allows one to compute the complete line profile (from the core to the far 
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Figure 26: Absorption cross section of the Ca I 4226 A line perturbed by neutral He atoms. 
The Lorentzian profile is compared to profiles obtained with the unified line broadening 
theory for different He densities relevant to the modeling of cool DZ white dwarfs [258). 


wings) of an atom interacting with other atoms or molecules in its environ- 
ment, and it remains valid up to the densities reached in the line-forming 
regions of most cool He-atmosphere white dwarfs. This theory is described 
at length in Refs. [256, 257]. It requires two important inputs: a potential 
energy curve and a transition dipole moment curve of the perturbing atom 
or molecule interacting with the radiating atom in its relevant states (i.e., 
the lower and upper states involved in the atomic transition being modeled). 

The unified line broadening theory has been very successful at modeling 
absorption lines in cool DZ white dwarfs [259, 260, 261, 262, 258, 37, 263], 
and its use is crucial to reach sound conclusions on the composition of rocky 
bodies accreted by those stars [264, 265]. Despite its success and impor- 
tance, this theory has so far only been applied to a small fraction of all the 
atomic transitions observed in cool white dwarfs. Most transitions are still 
modeled using Lorentzian profiles even in a regime where this is a very poor 
approximation (Figure 26). The chief reason for this state of affairs is the 
lack of appropriate potential energy and transition dipole moment calcula- 
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tions, which are obtained using high-accuracy quantum chemistry methods. 
Potential energy and transition moment curves for more transitions observed 
in cool white dwarfs would solve this problem and allow for a more robust 
determination of the composition of the very old planetary debris accreted by 
those stars. All C, Fe, Ni, Ti, Cr, K, and Li transitions observed in cool white 
dwarfs at visible and UV wavelengths are currently modeled with Lorentzian 
profiles and are therefore all good candidates for future work. 

One caveat of the unified line broadening theory is that it is based on po- 
tential energy and transition moment curves obtained in the infinite dilution 
limit (i.e., just one perturber interacting with the radiating atom). At the 
high densities characterizing the coolest DZ white dwarfs, this approxima- 
tion could become problematic as on average more than one perturber may 
be inside the interaction volume of the radiating atom and the interactions 
may not be linearly additive. This issue becomes particularly important 
when considering spectral lines produced by excited atomic states that are 
characterized by a larger spatial extent. 


5.8.8. Broadening by collisions with neutral particles: DB white dwarfs 

In the He-rich atmospheres of DB white dwarfs, broadening of neutral 
He lines by collisions with neutral He atoms becomes significant below Tig = 
16000 K. Even though the density remains low near the photosphere (ny. < 
10?! cm~?) and Lorentz line profiles are a good approximation to the line pro- 
file, the line width from collisions with neutral He atoms is poorly described 
in atmosphere models. This results in large uncertainties, increased scatter, 
and systematic deviations in the masses inferred from fitting DB white dwarf 
spectra for stars with Tyg ~ 11000 — 16000 K [266, 267]. The models use a 
combination of comparatively simple line broadening models [268, 269, 270] 
and ad hoc prescriptions. The widths based on the more advanced model 
of Ref. [271] are limited to T < 300K and not applicable to white dwarf 
atmospheres. Finally, the unified line broadening theory has been applied to 
broadening collisions with neutral He atoms but only for transitions in the 
UV and far red part of the spectrum [272] that are difficult to observe with 
ground based telescopes, if at all possible. There is a strong need for reliable 
calculations of the line widths of He broadened by collisions with He (and 
H) atoms for the transitions seen in the visible part of the spectrum. 
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5.4. Molecular opacities 


5.4.1. Cy absorption bands 

In cool white dwarf atmospheres, the formation of diatomic molecules can 
lead to observable molecular absorption bands in the optical spectrum. In 
general, those molecular bands are the result of a combination of electronic, 
vibrational, and rotational transitions. Bands from a few molecular species 
(e.g., CH, MgH) have been identified in white dwarfs but the most prominent 
features are by far the Cy Swan bands detected in cool DQ white dwarfs, 
which have atmospheres of He with traces of C. The Cy Swan bands are 
caused by the aI, 3 hs UB electronic transition and hundreds of different 
vibrational and rotational transitions. Figure 27 shows a typical cool DQ 
spectrum (in black), where the different vibrational bands are identified; the 
rotational transitions correspond to the finer structure within each of those 
vibrational bands. 

Current Cz line lists allow model atmospheres to precisely reproduce the 
Swan bands observed in most DQ white dwarfs [32, 273]. However, at very 
low temperatures (Tig < 6000 K), the C2 Swan bands are shifted and dis- 
torted compared to the predictions based on standard line lists [274, 275, 276]. 
Such objects are known as peculiar DQ, or DQpec, white dwarfs (an exam- 
ple is shown in gray in Figure 27). The cause of the DQ—>DQpec transition 
appears to be a density-driven shift of the electronic levels of the C2 molecule 
under the liquid-like densities that characterize DQpec white dwarfs. This 
interpretation is supported by DFT calculations that reveal a monotonic in- 
crease of the separation between the a®II,, and d®H, electronic levels with 
respect to the He density [277]. This results in a blue shift of the absorption 
bands. Once integrated over the densities of the atmospheric layers along 
our line of sight, this effect not only shifts the Swan bands, but also distorts 
them as indicated by the observations (Figure 27). 

While existing DFT calculations can qualitatively explain the DQpec phe- 
nomenon, the electronic energy transition shift is one order of magnitude too 
strong compared to what is needed to reproduce DQpec spectra |277, 273]. 
There are two possible explanations for this tension: (1) model atmospheres 
overestimate the atmospheric density of DQpec white dwarfs by an order of 
magnitude or (2) current DFT-based calculations of the Cy energy spectrum 
in dense He are inadequate. The first explanation appears unlikely since 
atmosphere models reproduce precisely the density-sensitive features of cool 
DZ white dwarfs in the same physical regime [192, 193, 32, 258, 262] and 
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Figure 27: Observed spectra of the DQ white dwarf BPM 27606 (black) showing the 
C2 Swan bands and of the DQpec white dwarf LP 315-42 (gray) where those bands are 
severely distorted. The spectra were taken from [31] and the Sloan Digital Sky Survey 
(SDSS) archive, respectively. 
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since the presence of H (which reduces the atmospheric density) is well con- 
strained by the absence or presence of a CH molecular band [32]. In practice, 
the distortion of Swan bands in DQpec models is included by simply linearly 
shifting the reference low-density Cy absorption spectrum. This approxima- 
tion is most likely too crude. For instance, the distortion of the electronic 
energy spectrum should in turn affect the vibrational energies. More work 
on the absorption of C2 in dense He is needed to solve this problem and to 
enable the reliable modeling of those peculiar objects. 

Another related question in need of further theoretical work is that of 
the effect of strong magnetic fields on Cy Swan bands. Fields of the order 
of 100 MG could provide another mechanism to distort Swan bands in white 
dwarfs [278], possibly explaining why some DQpec white dwarfs exist even 
in a regime where the photosphere is not particularly dense [273]. This effect 
remains poorly quantified. 


5.4.2. Collision-induced absorption 

In a dilute gas, the Hy molecule has no dipole moment and it produces 
only very weak quadrupole absorption features that are not detected in white 
dwarfs. At higher densities however, the Hj molecules undergo frequent col- 
lisions with other particles (Hz, H, and He in white dwarf atmospheres) that 
induce temporary dipole moments and lead to “collision-induced absorption” 
(CIA) [279]. In cool white dwarf atmospheres, this can give rise to broad and 
strong absorption features in the infrared [280, 134, 281]. With CIA being 
driven by collisions in the gas, it becomes stronger at increasing densities, 
and therefore it increases with decreasing effective temperature (Figure 16). 
In fact, in very cool white dwarfs, CIA is so strong at near infrared wave- 
lengths that it starts shifting the peak of the spectral energy distribution 
into the visible with decreasing temperature (at odds with the expectation 
from a simple blackbody spectrum) [195, 20]. White dwarfs that are past 
this turn-off point (i.e., when they start becoming bluer instead of redder 
as they age) are known as IR-faint (or ultracool) white dwarfs. These may 
be the oldest white dwarfs of our Galaxy, but their properties remain very 
uncertain [282, 283, 284, 20, 285, 286, 100]. 

Collision-induced absorption being a collisional process, its absorption 
coefficient can be expressed as a density expansion 


a(v) = a,(v)n + ao(v)n? + a3(v)n? +... (34) 


where n is the number density of the CIA absorber (we assume a single species 


83 


here for simplicity) and the a;(v) are the frequency-dependent absorption co- 
efficients involving 7 particles. By definition, a;(v) = 0 for a CIA process. 
Throughout the years, the accuracy of the theoretical CIA opacities used in 
model atmospheres has steadily increased [287, 288, 289]. However, those 
opacity calculations are normally based on a molecular scattering approach 
where only binary collisions (e.g., between Hj and He) and only ag(v) is 
considered. This is a valid approximation for cool H-dominated atmospheres 
(Figure 16), but it breaks down at the higher densities that characterize He- 
dominated atmospheres, where many-body collisions become important.'® 
The higher order coefficients are extremely difficult to evaluate with the tech- 
niques used to compute a2(v7). This problem has been addressed using DFT- 
MD simulations, and many-body collisions were indeed found to significantly 
alter the Hp-He CIA opacities for conditions relevant to cool He-atmosphere 
white dwarfs [290] (at least at \ = 1 um, CIA at shorter wavelengths has yet 
to be revisited). Similarly, it was shown that CIA can occur in cool He atmo- 
spheres even in the absence of Hy through three-body He-He-He collisions 
(a2(v) = 0 due to the symmetry of the He-He collision) [196]. Despite those 
advances, model atmospheres generally struggle to reproduce the observed 
spectral energy distributions of IR-faint white dwarfs (Figure 28). Recently, 
Ref. [286] presented model atmosphere calculations that markedly improved 
the agreement with observations. However, to achieve this result, they found 
that they had to resort to older CIA opacity tables [288]. Whether this is 
due to inadequacies in the current CIA opacities or to problems with another 
microphysical ingredient remains an open question. 


5.5. Phase diagrams of ultra-dense plasmas 


Modeling the fractionation (or phase separation) that takes place in crys- 
tallizing white dwarfs (Section 3.2) requires accurate phase diagrams of multi- 
component plasmas. Those phase diagrams yield the crystallization temper- 
ature of a given mixture and quantify the composition change between the 
coexisting solid and liquid phases upon crystallization. The problem of cal- 
culating those phase diagrams was first addressed in the 1980s and 1990s 
[121, 291, 123, 292], but here we will focus on three more modern techniques 
that have been used during the last decade. 


16Due to their high densities, even a trace of Hy in cool He-dominated atmospheres is 
sufficient to generate strong CIA features. 
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Figure 28: Best-fit model of the IR-faint white dwarf LHS 3250 [284]. The temperature, 
mass, and atmospheric composition of the model atmosphere were adjusted to obtain the 
best agreement possible with the observations (black error bars), but severe discrepancies 
remain. In most bandpasses, the discrepancies between the model and the observations are 
several times larger than the lo uncertainties. For non-IR-faint white dwarfs, a reduced 
chi-square 2 for such fits of ~ 1 is the norm; here y? = 75. Note that a blackbody 
spectrum at the best-fit temperature (Te¢ = 3060K) would peak at 1.7m. Here, the 
model spectrum peaks at much shorter wavelengths due to strong Hz CIA in the near- 
infrared. Note that Ref. [286] recently obtained better agreement with the observations, 
and found a much warmer Teg of 4990 K for this star. 
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5.5.1. Phase diagram calculation techniques 

A first phase diagram calculation method is the semi-analytic approach 
described in Ref. [293], which is based on analytic fits to the free energies 
of multi-ionic mixtures obtained with Monte Carlo simulations [294]. Ana- 
lytic fits allow a rapid determination of the phase diagram and makes this 
method applicable to mixtures with a large number of ionic species. Recent 
applications of this approach to white dwarf interiors include a calculation 
of the C-O-Ne phase diagram [125] and the study of the phase separation 
and settling of traces of actinides early in the crystallization phase [295]. Its 
main drawback is that it is only as accurate as the free energy fits on which 
it relies, and, unfortunately, those fits are unsatisfactory in several aspects. 
For instance, they do not include the effect of electron-ion screening (they 
are all based on the Monte Carlo or Molecular Dynamics simulations where 
the ions interact through bare Coulomb potentials); they implicitly assume 
that the vibration modes of the multi-component zero-temperature solid are 
identical to that of the one-component zero-temperature solid; and it is not 
always clear which functional form should be adopted for the free energy 
[294, 296]. Another issue with this approach is that those analytic fits are 
based on linear mixing rules that often require knowledge of the free energy 
of a strongly super-cooled liquid or super-heated solid (this problem is par- 
ticularly important for ionic mixtures with high charge ratios). Not only do 
current fits not extend into those regimes, but it is conceptually problematic 
to refer to a liquid one-component plasma at an extreme coupling regime 
that is well within the solid regime (e.g. I = 1000). This issue was recently 
discussed in Ref. [297]. 

A more direct approach that circumvents the need for analytic fits to the 
free energy is two-phase Molecular Dynamics simulations. Typically, those 
simulations are initialized with a solid and a liquid phase in contact with each 
other. The ions are allowed to diffuse through the solid—liquid interface, and 
eventually a state of equilibrium is reached. The compositions of the two co- 
existing phases then directly yield a point on the liquid-solid coexistence line. 
This approach has been successfully applied to mixtures relevant to white 
dwarf crystallization [75, 298, 299]. Unfortunately, the “first-principles” as- 
pect of this method comes with a very steep computational cost. Molecular 
Dynamics simulations that include a solid—liquid interface must be performed 
on very large systems to mitigate finite-size effects. As a consequence, only a 
very partial sampling of a given phase diagram can be achieved, which makes 
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this approach more suitable to verify the results produced by other, more ef- 
ficient methods than to calculate the detailed phase diagrams implemented 
in white dwarf codes. Applying this technique to study the fractionation of 
minor trace species such as ??Ne in a C-O white dwarf, is also challenging 
since a very large total number of particles is needed to include a reasonable 
number of ions of the trace species. 

Finally, a third approach is to directly integrate the Clapeyron relation 
that describes the shape of the liquid-solid phase boundary. When applying 
the Gibbs—Duhem relation 


sdT —vdP +> axidjy = 0, (35) 


where s and v are the entropy and volume per particle, respectively, and 2; 
is the number fraction of species 7, to the coexisting solid and liquid phases 
of a pure substance, we obtain the well-known Clausius—Clapeyron relation 


for the coexistence curve 
dP Ah 
(i) ER 36) 


where the derivative is along the coexistence curve and Ah and Av are the 
changes in specific enthalpy and specific volume across the transition. For 
mixtures, the Gibbs—Duhem relation gives rise to several forms of Clausius— 
Clapeyron equations. For a two-component system, we obtain a differential 
equation giving the slope dT'/djiz along the coexistence curve in terms of jug, 
Az, and Ah [300]. Given a starting point on the coexistence curve (usually 
2 = 0 or 1), the coexistence curve T(x2) can be obtained by integration at 
a fixed pressure. The (P,T,x2) phase diagram of a binary mixture can be 
generated by repeating the integration over a range of values for the pressure. 

This method has been used for some time in other domains [301, 300], 
but was only recently adapted to electron-ion plasmas [302]. Distinct Monte 
Carlo simulations for the two coexisting phases are used to evaluate the quan- 
tities needed to perform this integration. In the implementation of Ref. [302], 
all calculations are performed at constant pressure (the appropriate choice for 
phase transitions in white dwarfs) and explicitly include the electron back- 
ground (which is usually neglected with other methods where the constant- 
volume approximation is ordinarily made). Just like the Molecular Dynamics 
approach, the Clapeyron method avoids the need for thorny analytic fits of 
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Figure 29: Carbon-Oxygen phase diagram. The solid lines are the liquidus (above which 
the plasma is liquid), and the dashed lines are the solidus (below which the plasma forms 
a stable solid). The temperature is given in units of the melting temperature of a pure 
C plasma, and the horizontal axis gives the O number concentration (pure C at ro = 0 
and pure O at x5 = 1). The separation between the solidus and the liquidus at constant 
temperature gives the composition change upon solidification of the plasma. Note that 
this phase diagram is virtually independent of the pressure under white dwarf conditions 
because the pressure is entirely dominated by the degenerate electron gas and therefore 
decoupled from the ions. The results showed here are taken from Refs. [75, 293, 126]. 


the free energy. But it is also much less computationally costly than the two- 
phase Molecular Dynamics technique. This is mainly due to the fact that 
the two phases are simulated separately. This means that there is no solid— 
liquid interface in the Monte Carlo simulations, thereby greatly reducing the 
number of particles needed to mitigate finite-size effects (from 10,000s of ions 
to just under 1000). This approach has therefore allowed the calculation of 
highly accurate and well-sampled phase diagrams for a variety of two- and 
three-component mixtures relevant to white dwarfs interiors [126, 17, 303]. 
Figure 29 shows three versions of the C-O phase diagram, obtained using 
the three techniques just described. It is reassuring to see that these three 
totally independent approaches lead to very similar results. In particular, all 
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three methods predict an azeotropic phase diagram with the azeotropic point 
at %o & 0.2 (inset). Note that for this mixture, the semi-analytic approach 
is not expected to be affected by the super-cooled liquid (and super-heated 
solid) problem, due to the small charge ratio (Zo/Zc = 4/3) of the mixture. 
The coarseness of the Molecular Dynamics-based phase diagram makes it 
inappropriate for direct use in white dwarf codes, but is nevertheless useful 
to verify the results of the other methods. 


5.5.2. Some recent results: Crystallization of trace species 

The excellent agreement shown above for the two-component C-O phase 
diagram suggests that this problem is mostly solved. Accordingly, the most 
recent works on phase diagrams of white dwarf interiors have instead focused 
on the role played by trace species. Even if C and O largely dominate most 
white dwarf cores, the fractionation of traces species can play a dispropor- 
tionate role on the evolution of the star. Three recent works demonstrate this 
point and show that there are still important physical processes to uncover 
by calculating phase diagrams for mixtures found in white dwarfs interiors. 

Ne distillation — Ref. [17] investigated the fractionation of ??Ne in 
C-O white dwarfs !” and showed that the solid phase can be significantly 
depleted in ??Ne compared to the coexisting liquid. Because 7*Ne has a 
higher mass per electron than !°O and !2C, the ?*Ne-poor crystals are lighter 
than the liquid in which they form. This is expected to lead to a “distillation” 
process [122] that can transport large quantities of ??Ne toward the central 
layers. This in turn liberates the gravitational energy stored in ?*Ne and can 
lead to multi-Gyr cooling delays that appear to be needed to explain recent 
astronomical observations [120, 13, 20]. 

Fe-rich crystallite formation — °°Fe is another neutron-rich trace 
species present in white dwarf cores (with a mass fraction of the order of 
0.1%). Ref. [127] found that °°Fe should separate into Fe-rich crystals be- 
fore the rest of the ionic mixture freezes. Those heavier crystals should sink 
rapidly to the center of the white dwarf, thereby releasing gravitational en- 
ergy and slowing down the cooling of the white dwarf by a few hundred 
million years [304]. 

Actinide fractionation — Ref. [295] looked at the case of even lesser 
species, actinides. Because of their much higher ion charge, resulting in 


17 22Ne makes up 1 — 2% of the mass of C-O cores. 
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strong plasma coupling (Equation 8), actinides crystallize at temperatures 
well above that of C and O. They argue that when the white dwarf begins 
to crystallize, the solid phase could be so enriched in actinides that it could 
support a fission chain reaction, thereby potentially igniting C burning and 
leading to a supernova explosion. While this scenario needs to be supported 
by more calculations, the prospect of a completely new supernova mechanism 
is tantalizing. 


5.5.8. Remaining challenges 

The three methods described above treat the ions as classical particles. 
This approximation is valid in the limit where the mean interparticle dis- 
tance is larger than the thermal de Broglie wavelength of the ions, but it 
starts to break down in white dwarf cores (Figures 10a and 12b) [305, 306]. 
Future work should therefore investigate how quantum effects impact phase 
diagrams under white dwarf conditions. Very recently, Ref. [297] included 
the quantum behavior of ions in their equation of state and find a very weak 
effect on the crystallization temperature, with [’ decreasing from 175 to 168 
for carbon. However, they could not study the impact of quantum effects on 
fractionation due to limitations in their EOS implementation. 

Another important approximation that is made by all three methods 
listed in Section 5.5.1 concerns the lattice structure of the solid phase. For 
a one-component plasma, it is well established that the most stable solid 
structure is a bcc lattice [307]. However, in multi-component plasmas, other 
more stable configurations may be possible [294, 308, 309, 310]. It remains 
unclear to what extent the bcc lattice is a good approximation of the real 
crystalline structure, especially for highly asymmetric mixtures. Moreover, 
the solid may well not even form a regular lattice structure but could instead 
form a more complex polycrystalline mixture. 

Finally, reliably modeling the fractionation of high-Z trace species re- 
mains challenging. Molecular Dynamics and Monte Carlo-based methods 
demand simulations with a very high total number of ions in order to include 
a reasonable number of ions of the high-Z trace species, and the semi-analytic 
approach is precarious due to the limitations of the analytic free energy fits. 


5.6. Diffusion coefficients for trace species 


In most instances of diffusion in white dwarfs, we are concerned with 
the inter-diffusion in a binary mixture of ions. For a mixture with number 
concentrations 7; and x2 = 1— 4, arbitrary electron degeneracy and plasma 
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coupling, the equation for the relative velocity of the ion species wiz due to 
diffusion is 


D ce ) dln x (A= AN) tos ( Ly — Zy joe 
Wi = = 
cael aaa dr Fea GI NE Aaa) aie 


| 


Ar 


dr 

(37) 
where Dj» is the coefficient of inter-diffusion, A; and Z; the atomic mass 
(in a.m.u.) and charge of ions of species i, P; the total ionic pressure, mo 
the atomic mass unit, y = X2/x,, and a, is the (dimensionless) thermal 
diffusion factor [6, 144]. The first term on the right-hand side of Equation 37 
corresponds to “ordinary” chemical diffusion. For trace elements (12 < 21), 
this term is negligible in white dwarfs [311]. The second and third terms 
describe barodiffusion, known as “gravitational settling” in stars, caused by 
the pressure gradient associated with the star’s gravitational field and the 
induced electric field that maintains local charge neutrality for a mixture 
of charged species of different masses in a gravity field. This is generally 
the dominant term in white dwarf envelopes [142]. The last term is the 
contribution of thermal diffusion. The coefficient of thermal diffusion D, = 
Dj,2ar is determined by ion-ion and electron-ion collisions and can be written 
aS Ap = AyQ+ Ae + Q2-. The contributions of aj. and a2. can be comparable 
or even much larger than the ion-ion term aj. [142]. 

When Equation 37 is applied to a fully ionized plasma, the value of Z; is 
unambiguous. However, it is common for one of the elements to be partially 
ionized [144], in which case a separate equation is necessary for each charge 
state, treated as a separate species. An average diffusion velocity for a given 
element can be defined for such a multi-component diffusion problem but it 
is a common approximation to use a single diffusion equation for an element 
with an average ion charge instead [312, 145, 144]. 


5.6.1. Diffusion in the envelope 

In stellar astrophysics, the evaluation of diffusion coefficients is usually 
based on the expression derived from the solution of the Boltzmann kinetic 
equation [313, 314]. The resulting transport coefficients are expressed in 
terms of collision integrals involving scattering cross sections for isolated, 
binary collisions. In systems where particles interact with a pure Coulomb 
potential, these integrals diverge. In a real plasma however, the pair interac- 
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tions are screened by correlations between particles, leading to a short-range 
effective pair potential and finite collision integrals. An effective pair poten- 
tial is appropriate if the dynamics of charged particles can be described as 
a succession of binary collisions in this effective potential, a condition that 
corresponds to a weakly coupled plasma (I. < 1). The simplest form of ef- 
fective interaction is a pure Coulomb potential with an upper radius cutoff, 
usually chosen to be the Debye screening length Ap. A more physically cor- 
rect description of the pair potential uses the statically screened Coulomb 
(or Yukawa) potential [315] 

Viel) = Berl , (38) 
which is found to be applicable for [ < 0.3. The most common source of 
collision integrals used in white dwarf models are the fits of Ref. [143] that 
use the Yukawa effective potential with a screening length given by 


A= Wax) hy Gl (39) 


where Ap, is the Debye screening length and a is the ion sphere radius. This 
choice provides a heuristic correction to the screening length that extends 
the applicability of the approach into the strong coupling regime. Ref. [316] 
improved on this calculation by accounting for electron degeneracy in Ap and 
substituting a smooth, physically motivated interpolation form for Equation 
39. Note that the models of Refs. [143] and [316] require a separate input of 
the ion charges Z; and Z. In white dwarf envelopes, the dominant ions (H 
and He) are usually fully ionized but trace species are typically low-Z ions 
(Z ~ 6 — 26) that are partially ionized. The calculation of the charge state(s) 
of such ions in the warm dense matter regime is an ongoing challenge [317]. 

More recently, a rigorous extension of the Boltzmann kinetic equations 
into the moderate coupling regime, called the effective potential theory has 
been developed. It can be shown that the proper effective potential is the 
ion-ion potential of mean force Vyy(r), defined as 


g(r) = exp(—Vur(r)/keT) , (40) 


where g(r) is the pair distribution function of the ions [318]. The potential 
of mean force includes the direct binary interaction between a pair of parti- 
cles V(r) and corrections that account for the presence of the surrounding 
correlated particles. This formalism is general and can be applied to the 
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calculation of collision integrals for a system with any direct pair potential 
(e.g. Equation 38), as long as g(r) can be calculated or experimentally mea- 
sured. The effective potential theory has been extended to plasma mixtures 
[319, 320] and applied to thermal diffusion [321]. It accounts for ion screening 
through V\;7(r) while the electron screening is implicit in the direct pair po- 
tential V(r), as the electron contribution to the screening length in Equation 
38 for example. The principal remaining difficulty is to develop a realistic 
model for electron screening in regimes where the Yukawa potential is not 
appropriate. 

Average atom models explicitly solve for the quantum mechanical states 
(bound and free) of the electrons surrounding a nucleus, including the effects 
of the surrounding plasma. While approximate, average atom models nat- 
urally account for pressure ionization and strongly non-linear screening (see 
[322] for a recent overview). They provide self-consistent ion charge, elec- 
tron cloud spatial distribution and the ion-ion interaction potential. They 
have been very useful for the description of states of warm dense matter, 
including the EOS and transport properties. Recently, Ref. [323] combined 
the potential of mean force extracted from an average atom model [324] with 
the effective potential theory to compute interdiffusion coefficients and ionic 
thermal diffusion coefficients (a@j2) of Si and Ca under the conditions found 
in white dwarfs. In the weakly coupled regime, which is relevant to all H-rich 
white dwarfs and those He-rich white dwarfs hotter than Tig ~ 14000 K, the 
results agree with those of Ref. [143] to within 25%. In the cooler He-rich 
stars, however, the background He plasma becomes strongly coupled and the 
differences grow to a factor of 2.4, which is significant (Figure 30). These 
values are obtained when the ion charges calculated with the average atom 
model are used as input to the Ref. [143] formalism. With the charge state 
obtained with a more heuristic ionization model [325], the coefficients of in- 
terdiffusion of Ref. [142] can change by as much as 30% across all coupling 
regimes. This shows that a realistic ionization model to compute diffusion 
coefficients is as important as using an advanced theory of transport coeffi- 
cients [144, 326]. This effect is compounded when the ratio of interdiffusion 
coefficients of different elements is the more relevant quantity, such as in the 
estimation of the elemental composition of the accreted material onto pol- 
luted DZ white dwarfs. It is also important when computing the upward 
diffusion tail of C at the C-O/He interface when modeling the dredge up of 
C observed in DQ and DQpec white dwarfs. 

The state of the art in the calculation of transport properties are ab ini- 
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Figure 30: Interdiffusion coefficient of a trace of Ca in He for the conditions at the bottom 
of the convection zone of He-rich white dwarf models as a function of the model Tag. The 
(p,T) conditions vary as a function of T.g. The values from the Effective Potential Theory 
combined with the average atom model of the plasma are shown in blue (“EPT”). The 
black symbols (“PAMD”) show the values obtained with classical Molecular Dynamics 
simulations using the direct potential from the same average atom model. Ideally, these 
would agree perfectly with the EPT result. The differences and scatter in these data show 
the difficulty of obtaining converged molecular dynamics results for the diffusion of a trace 
species. The simpler model of Paquette et al. [143] with the same ion charges as the 
EPT and PAMD calculations is shown with the dotted red line. Changing the ion charges 
to those from a simpler, somewhat heuristic ionization model [325] results in the solid 
red line. In this case, it fortuitously improves the agreement with the EPT calculation. 
Adapted from Figure 5 and Table 3 of Ref. [323]. 
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tio computer simulations of plasmas that consider only nuclei and quantum 
mechanical electrons. These methods involve the fewest number of approxi- 
mations but are computationally very expensive. It is difficult to obtain con- 
verged values for interdiffusion coefficients when considering a trace species 
or if the mass ratio or charge ratio between the two species is large. At 
present, ab initio simulations are impractical to generate the large tables 
of diffusion coefficients for the many different elements of interest to model 
white dwarfs. 

To summarize, the combination of an average atom model with the Ef- 
fective Potential Theory [323] is a significant advance in the calculation of 
interdiffusion coefficients compared to the approach of Ref. [143] and vari- 
ations thereof (e.g. Ref. [316]). These improvements matter primarily for 
the cooler He-rich white dwarfs polluted with low-Z elements (type DZ) and 
those with traces of dredged-up C (spectral types DQ and DQpec). The 
contribution of ion-electron interactions to the coefficient of thermal diffu- 
sion (@j- and a2.) remains to be evaluated with a comparable level of theory 
for quantum mechanical electrons (e.g. Ref. [327]). Benchmark ab initio 
calculations (e.g. with Density Functional Theory - Molecular Dynamics 
simulations) for a few representative cases would be very valuable to validate 
the more efficient but approximate models that can be applied realistically 
to the extensive computations of diffusion coefficients. 


5.6.2. Diffusion in the core 

Diffusion (or gravitational settling) of trace species in the core is also an 
important process. Neutron-rich species (with A > 2Z) experience a net 
downward force in white dwarf cores and therefore tend to diffuse toward the 
centre (wi2 < 0 in Equation 37). In C-O white dwarfs, the most important 
species to undergo this process is 77Ne, which makes up ~ 1 — 2% of the 
core mass. This transport releases a sizeable amount of gravitational en- 
ergy, which slows down white dwarf cooling [328, 329, 330, 331]. This effect 
must be considered in the building of white dwarf cooling sequences used for 
age-dating applications, and its efficiency depends on the relevant diffusion 
coefficients. 

Gravitational settling of ??Ne in the core takes place in a fully ionized 
plasma, eliminating some of the complications discussed above for the dif- 
fusion of low-Z elements in the envelope. Many works have focused on 
determining self-diffusion coefficients in a one-component Yukawa plasma 
[332, 333, 334, 335, 336, 337] and this quantity is now well known. How- 
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ever, the inter-diffusion coefficients needed to model diffusion in the multi- 
component plasma that makes up the core are less well constrained. Linear 
mixing rules are often used, but their accuracy should be further validated. 


5.7. Thermal conductivity 


The calculation of the electron thermal conductivity in stellar astrophysics 
and in white dwarfs in particular has a long history, most recently sum- 
marized in Ref. [338]. The current standard source of conductive opacity 
used in white dwarf models is Ref. [138]. In the C-O core of white dwarfs, 
conductive opacities accounting for the strong electron degeneracy, strong 
plasma coupling and partial relativity (Figures 9, 10, 15) are well under- 
stood. The regime of partial degeneracy and moderate coupling found in the 
He and H layers of white dwarfs is more challenging. The simple Spitzer— 
Harm expression [339] for the electron conductivity in a weakly coupled, non- 
degenerate and non-relativistic plasma can be extended to stronger coupling 
with the Chapman—Enskog solution of the kinetic equation for transport us- 
ing the Landau-Fokker—Planck collision operator. The recent development 
of a quantum mechanical (accounting for electron degeneracy) form of the 
Landau-Fokker—Planck kinetic equation for plasmas (qLFP; [340]), combined 
with the potential of mean force from an average atom model [341] to de- 
scribe electron-ion and electron-electron collisions, has allowed the extension 
of calculations of the electron thermal conductivity into the partially degener- 
ate regime where both electron-electron and electron-ion collisions contribute 
to a comparable degree [342]. Previously, the electron-electron contribution 
was interpolated in this regime [138]. The conductive opacities for H and He 
in white dwarfs computed with this new model [139] depart from those of 
Ref. [138] in a narrow region of intermediate degeneracy with 0 ~ 0.1 - 1 
(Figure 31). The lower conductive opacities result in faster cooling of white 
dwarf models. A typical 0.6 Mo H-rich white dwarf reaches the end of the 
cooling sequence (Tig * 4000 K) about 1 Gyr earlier, while a massive 0.9 Ms 
H-rich white dwarf will reach the same temperature about 2 Gyr sooner. The 
correction in cooling time is about 3 times smaller for He-rich white dwarfs, 
owing to the smaller change in the conductive opacity of He. Given that the 
cooling age of such stars is 8-10 Gyr, these corrections are substantial. The 
correction to the thermal conductivity due to electron-electron collisions in 
higher Z plasmas such as carbon will be smaller than for He and are likely 
negligible. 
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Figure 31: Contours of the ratio of the standard conductive opacity «1° [138] to the more 
recent opacities K¢4FP [139]. The long axis of the green band corresponds to 0 ~ 0.3 for 
both H (left panel) and He (right panel). Adapted from Ref. [139]. 


While Ref. [139] limits their calculation to a range of validity estimated 
to be O = 0.1, Ref. [338] points out that the bulk of the change in the 
conductive opacity calculated by Ref. [139] occurs for © < 1 where the qLFP 
model is less accurate and likely departs too strongly from the standard 
conductive opacities of Ref. [138]. Exploring a range of possible behaviors 
of the conductive opacity for 0 < 1 while keeping the opacities of Ref. [139] 
for 0 > 1, they find that the effect on the cooling time of old white dwarfs 
can range from slightly shortened to close to the large values reported by 
Ref. [139]. This recent work points to the need for a better model for the 
calculation of the thermal conductivity of H and He for 0 ~ 0.1 —1. 


5.8. Laboratory experiments 


The extreme physical conditions encountered in white dwarfs are very 
challenging—when not impossible—to reproduce in the laboratory. As a re- 
sult, there is an extreme paucity of experimental data that is directly relevant 
to white dwarfs. Much of the constitutive physics that is input to the models 
is based on physical theories and models that are of unknown accuracy’. 


‘8Tnterestingly, the problem of the lack of experiments under white dwarf conditions 
can sometimes be inverted by using white dwarfs as laboratories to test physical theory 
in otherwise inaccessible regimes such as the crystallization of the core (Section 5.5) [343, 


ot 


Contemporary white dwarf research is increasingly concerned with their de- 
tailed properties and the stellar models are confronted with a very rich trove 
of observations that calls for increased physical realism and detail. Much 
of white dwarf astrophysics relies on our confidence in the validity of the 
underlying physical models. 

Unique, world-class experimental facilities built to create and study mat- 
ter at extreme conditions, with a focus on nuclear fusion research, have be- 
come available for research of a more academic nature.!? Ongoing exper- 
iments at these facilities study various properties of white dwarf plasmas. 
This new, exciting development in white dwarf science heralds a fresh eval- 
uation of some of the constitutive physics and provides opportunities for 
original experimental concepts. 


5.8.1. Line profiles 

Experiments that study the Stark effect in H plasmas can achieve tem- 
peratures and electron densities found in white dwarf atmospheres (T ~ 1 — 
3eV, ne ~ 5 x 101° — 10?°cm~3) with a variety of techniques such as Laser 
Induced Breakdown Spectroscopy, gas-liner Z-pinch, and arc discharge. The 
vast majority of these studies are limited to measurements of a single emis- 
sion line (typically Ha or H’) and often focus on the Stark line width and 
line shift [344, 345, 346, 347]. On the other hand, accurate modeling of the 
absorption lines in white dwarf spectra requires detailed theoretical line pro- 
files for several lines in a series up to the merging of the pressure-broadened 
lines near the limit of the series. For five decades, the main experimental 
reference for Stark broadening of the Balmer lines of hydrogen under condi- 
tions similar to those in white dwarf atmospheres has been Ref. [348], who 
measured the line profiles of Ha to He in emission as well as the merging of 
lines near the Balmer limit. 

New line profile experiments are underway at the Z-pinch pulsed-power 
facility at Sandia National Laboratories. The Z machine is a unique and 
versatile experimental facility that can create extreme physical conditions 
for a wide variety of high energy density experiments including equation of 
state, opacities, and magnetic confinement fusion [349]. This facility consists 


126, 17]. 

19Two programs that provide access to such facilities to the broader community are the Z 
Fundamental Science Program at Sandia National Laboratories and the National Ignition 
Facility Discovery Science program at the Lawrence Livermore National Laboratory. 
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of a bank of capacitors that can deliver an electrical pulse of up to 30 MA 
on time scales of ~ 100 ns onto cm-scale targets. It can be used to generate 
a hot plasma emitting a 330 TW flash of soft X-rays. A new experimental 
platform was developed to use this X-ray source to heat and ionize hydro- 
gen in a large (~ 10cm) gas cell. Conditions achieved in the H plasma are 
ne ~ 0.3 -9 x10! cm? and T ~ 1 - 1.5eV which are typical of DA white 
dwarf atmospheres (Figure 32). Simultaneous spectroscopic measurements 
of the HG, Hy and Ho lines in both emission and absorption are taken with 
a streak camera which allows tracking of the time evolution of the plasma 
conditions [350, 351]. The experiment is designed to minimize spatial inho- 
mogeneities in the plasma. Separate fits of the HG and the Hy lines observed 
in absorption give inconsistent values of n., however [351], which echoes the 
results of fitting different Balmer lines in white dwarf spectra. These results 
have prompted refinements in the experimental design and a new effort in 
modeling line profiles in low-Z elements (Section 5.3.1) [248, 250, 249, 352]. 
This experimental approach has now been extended to the neutral He lines 
5015A and 5875 A [353] that are among those used in characterizing DB 
white dwarfs. Because of the large difference in the ionization potentials of 
H and He, it is possible to probe both the broadening of He lines by colli- 
sions with neutral atoms and the Stark broadening caused by the surrounding 
plasma by mixing H (a source of electrons and protons) and He (a source 
of neutral atoms) by controlling the experimental conditions and the H/He 
mixing ratio. This is a particularly valuable approach in view of the uncer- 
tainties in the neutral broadening of He lines (Section 5.3.3). Measurements 
of C and O lines relevant to hot DQ white dwarfs are also under way [354]. 


5.8.2. Equation of state and K-shell conization 

A very recent development is the capability to produce plasmas at con- 
ditions found well below the atmosphere of white dwarfs. The National Ig- 
nition Facility (NIF) is the most energetic laser facility in the world and was 
designed to achieve nuclear fusion by inertial confinement of a deuterium- 
tritium mixture. The 192 laser beams deliver up to 1.8MJ of energy on a 
millimeter-size target on a time scale of a few nanoseconds [356]. Typically, 
the UV lasers shine on the inside surface of a small cylindrical hohlraum 
that produces thermal X-ray radiation that fills the cavity. The target is a 
sphere of ~ 1mm diameter suspended at the center of the hohlraum. Its 
outer layers are ablated by the X-ray radiation, sending a strong converging 
shock wave towards its center which rapidly compresses to a few 100 Gbar of 
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Figure 32: Electron density along pure hydrogen white dwarf atmosphere models of Tog 
= 4000 K to 20000K and log g = 8 (from left to right). The depth in the atmosphere is 
represented by the temperature, which increases downward in the star. The solid dot on 
each curve indicates the location of the photosphere, the nominal level where the spectrum 
is formed. The purple and blue rectangles represent the experimental conditions achieved 
by Ref. [348] and at Z [350, 355]. 
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Figure 33: Sequence of temperature-density structures of a cooling hot DQ white dwarf 
of 1 Mo from Teg = 24000 (bottom profile) to 18000K (top profile). Upon cooling, 
the structure evolves from right to left. Regions of convective mixing are shown in red. 
Calculated Hugoniots of diamond (dark blue) and oxygen (light blue) are shown [358]. 
The conditions to be achieved in the NIF Gbar experiments are highlighted by the green 
circle, and the conditions reached by Ref. [359] are represented by the box on the left. 


pressure and heats up to several keV. Inertial confinement fusion (ICF) tar- 
gets contain a mixture of deuterium and tritium, and fusion reactions occur 
in the compressed target [357]. 

This basic concept of spherical shock compression within a hohlraum has 
been adapted to conduct equation of state measurements at pressures ap- 
proaching a billion atmospheres (Gbar). The Gbar platform on NIF uses 
X-ray streak radiography to monitor the time evolution of the spherically 
imploding capsule, the convergence of the shock front and the density con- 
trast between the unshocked and shocked material. Combined with the 
Rankine-Hugoniot shock conditions, these measurements provide a time- 
resolved record of the pressure and density along the principal Hugoniot well 
into the warm dense matter regime. The Gbar platform has produced precise 
Hugoniot data that can discriminate between different equation of state mod- 
els at unprecedented pressures up to 400 Mbar for CoHi9 [360] and 800 Mbar 
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for diamond. Pressures of ~ 1 Gbar on the principal Hugoniot correspond to 
the ionization of the K electronic shells of C and O. This is a regime that is 
critical to many white dwarf processes (convection and mixing, the driving 
of pulsations, diffusion) and most challenging to model. Yet, the pressure 
in EOS tables used in white dwarf models can differ by 40% in this regime. 
State-of-the-art equations of state [358] remain untested under those condi- 
tions. In a ongoing experimental campaign, the Gbar platform on NIF is 
used to measure the shock Hugoniot of C and O well into the 1 Gbar regime 
to constrain models of the partial ionization of the K shell electrons. Typical 
conditions will be T ~ 100 —200eV, P ~ 0.1—1.5 Gbar and p up to 5gcm~? 
(oxygen) and 15gcm~? (carbon). These conditions correspond to the base 
of the convection zone of hot DQ white dwarfs (Figure 33). The prospect of 
experimentally probing the conditions higher up in the envelopes of hot DQs 
(at lower pressures and temperatures) was recently explored by Ref. [359] 
using isochorically heated carbon foams. 

This new development opens interesting possibilities for other experi- 
ments at those extreme conditions. The charge of a low-Z trace element 
that is gravitationally settling in a fully ionized background of H or He is a 
key factor in determining its diffusivity [323]. Under the conditions of inter- 
est, elements such as Mg, Si, Ca, Fe are partially ionized and the equation of 
state models used to obtain their average charge make widely different pre- 
dictions. An experiment that could measure the charge state of ions under 
white dwarf or similar conditions to challenge these models would be very 
valuable. This could be achieved by measuring the free electron density n- 
if the background species (H or He) is fully ionized. While the average ion 
charge in a warm dense plasma is an ill-defined parameter [317], the relevant 
quantity for ionic inter-diffusion is the effective (screened) ion-ion interac- 
tion, which can be probed experimentally by measuring the pair distribution 
function or the structure factor. A direct measurement of diffusivity would 
be even more valuable but is likely far more challenging. 


5.8.8. Opacities 

We have seen that the opacity in cool, pure He atmospheres of spectral 
type DC presents difficulties from the point of view of the underlying chemical 
equilibrium, particularly in estimating the number density of free electrons 
(Section 5.2.4). Similarly, the cross-sections of all four important absorption 
processes in He atmospheres (He Rayleigh scattering, He~ free-free, Hey 
bound-free, and triple-He collision-induced opacity, Figure 18) are affected 
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by the correlations in fluid-density atmospheres [361, 196]. Opacity measure- 
ments are very much necessary to validate the models for both the chemical 
equilibrium and the frequency dependence of these continuum absorption 
cross-sections. This is particularly important in DC white dwarfs because of 
the absence of spectral lines that provide powerful diagnostics of conditions 
in the atmosphere. Of particular interest are laser-heated diamond anvil cell 
experiments that measured the transmittance of heated and compressed He 
[362]. The data were obtained along two isobars (22 and 52GPa) at tem- 
peratures ranging from ~ 3200 to 16000 K,”? which match the conditions in 
atmospheres of Tg = 4000 and 5000K (Figure 16). The data show a rapid 
transition from a transparent insulating fluid to an opaque poor conductor 
around 10000 K — a lower temperature than expected from models. In addi- 
tion, the absorption coefficient increases toward shorter wavelengths, which 
is opposite to the trend of the dilute gas He~ free-free cross-section (Figure 
18). These are important results that would affect the predicted spectral 
energy distribution of He-atmosphere DC white dwarfs. Unfortunately, the 
range of data is too narrow in wavelength and pressure to confidently develop 
and calibrate a model that can be implemented in white dwarf atmosphere 
models. Additional He data with this technique would be very helpful. 

A well-known problem in stellar astrophysics is the tension between the 
increasingly reliable determination of the composition of the Sun that leads 
to interior models that do not agree with the exquisite data of helioseismol- 
ogy [363]. An increase in the opacity of the solar plasma could resolve this 
discrepancy. A considerable experimental effort has been invested in measur- 
ing the opacity of iron (the principal suspect) at the conditions found at the 
radiative/convective boundary in the Sun [364]. The experiments conducted 
at the Z machine give surprising results that have motivated a very detailed 
series of experiments and theoretical work in atomic physics [3865]. They have 
also prompted the development of the new Opacity platform at the NIF to 
investigate the opacity of Fe with a completely different technique [366]. This 
Opacity platform is also being applied to white dwarf physics to measure the 
opacity of O in the partial ionization regime, reaching temperatures of 100— 
200eV and densities of 0.01-1gcm~* (Figure 13). The opacity is difficult 
to model in this dense, partially ionized regime. These X-ray opacity mea- 
surements, spanning photon energies of 0.6 to 1.6keV, will provide strong 


20The corresponding densities are ~ 0.5 to 1.2g¢cm7°. 
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constrains on the theory. 

The Z pulsed-power machine and the NIF are unique, world-class facilities 
that complement powerful “table top” techniques like laser-heated diamond 
anvil cells. Together, they allow for the first time direct measurements of 
equation of state, opacities and line profiles of key materials (H, He, C and O) 
under the conditions found in white dwarf stars and bring strong constraints 
on untested theoretical calculations. The challenging experiments described 
above are the first efforts in this new area of laboratory astrophysics. We 
anticipate that during the next decade, more advanced techniques and cre- 
ative approaches will be developed at these and other facilities with powerful 
lasers [367, 368, 369], X-ray free electron lasers [370, 371, 372], and the next 
generation of pulsed-power machines [349] to probe white dwarf matter in 
new ways. 


6. Conclusion 


Astrophysics is a vast playground for the application of physical theory in 
exotic systems and challenging regimes. Stars of all types have been studied 
for well over a century and we understand them well enough to ask very 
detailed questions whose answers require close attention to the constitutive 
physics of the models. This is particularly true of white dwarfs where the 
physical conditions are quite different from those in other types of stars. 

As we have seen, there has been steady progress in the modeling of the 
physics of white dwarfs—particularly over the past decade and a half—that 
has resulted in tangible improvements in white dwarfs models, better agree- 
ment with the observations, the clarification of astrophysical puzzles and, 
in a few instances, their solution. Yet, many problems remain untouched 
or unsolved, or the advances, while tangible, are not sufficient to match the 
observations. In other instances, the physical theory has been developed but 
not yet applied to white dwarf modeling. It is important to keep in mind 
that state-of-the-art stellar models are the product of complex multi-physics 
codes and that the adoption of a better physics model requires that a prac- 
tical implementation be possible. 

As a summary and for convenience, we list the various current problems 
that we have discussed with references to the relevant sections. 


e Application and experimental validation of Zeeman—Stark line broad- 
ening and line shifts in the non-perturbative regime to magnetic white 
dwarfs (Section 3.6) 
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Ionization and opacity of He at p ~ 0.1—1gcem™? and T~ 0.1—1eV 
(Section 5.2.4) 


Dissociation equilibrium of diatomic molecules (Hy, C2, He}, HeH*) in 
dense He (Section 5.2.3) 


Line broadening of all the most important transitions in cool metal- 
polluted white dwarfs (Section 5.3.2) 


Broadening of optical lines of He by collisions with He atoms (Section 
Dae) 


Modeling and measurement of the distortion of Cy molecular bands in 
dense He and in the presence of magnetic fields (Section 5.4.1) 


Ion quantum effects on phase diagram calculations (Section 5.5.3) 


Crystallization dynamics and lattice structure in crystallized white 
dwarfs for different mixtures of species including impurities (Section 5.5.3) 


Robust calculations of high-Z (such as uranium) impurity fractionation 
in crystallizing white dwarfs (Section 5.5.3) 


Inter-diffusion coefficients: Benchmark ab initio calculations and new 
tables of Diz and a; based on modern theory, experimental measure- 
ment of the charge state of low-Z impurities in a H or He plasma (Sec- 
tion 5.6.1) 


Calculations of the thermal conductivity of H and He in the 0 ~ 0.1 — 
1 regime (Section 5.7) 


Experimental measurements of the absorption coefficient of dense, heated 
He under the conditions in cool He-rich DC atmospheres (Section 5.8.3) 


Experimental validation of the latest Stark line shapes for H and He 
(Section 5.3.1) 


Experimental measurement of the H2-He collision-induced absorption 
in cool white dwarf atmosphere conditions (Section 5.4.2) 


Line formation in the context of an inhomogeneous convective atmo- 
sphere (Section 5.1) 


105 


e 3D hydrodynamics simulations of convective overshoot and thermoha- 
line instabilities in the context of chemical mixing (Section 5.1.1) 


e Magnetic effects on convection and envelope structure, as well as verti- 
cal diffusion and horizontal spreading of accreted material (Sections 5.1.2 
and 5.1.3) 


This long enumeration is by no means exhaustive. Many other prob- 
lems arise in other types of white dwarfs that are not discussed herein. For 
instance, non-local thermodynamic equilibrium (NLTE) prevails in the at- 
mospheres of hot white dwarfs with Tig 2 40000 K. The characterization of 
those stars, particularly their composition, is hampered by our incomplete 
knowledge of the relevant atomic physics. In young white dwarfs, stellar 
winds driven by radiation gently blow out the surface layers, compete with 
diffusion and alter the surface composition. Modeling of these winds also 
depends on atomic physics data, some of it poorly known. White dwarfs 
are commonly found in binary star systems where matter from the other 
star is transferred onto the white dwarf. This leads to many phenomena not 
seen in single white dwarfs. The transferred material may trigger nuclear 
fusion under a wide variety of scenarios. If the accreting white dwarf has a 
strong magnetic field, the infalling matter is funneled to the magnetic poles 
in two accreting columns. Because of their very high density, white dwarfs 
have been proposed as potential “detectors” of various dark matter particles 
whose accumulation in their core could affect the cooling rate. The breadth 
of interesting physics associated with white dwarfs is quite remarkable. 

The challenges that we have explored in some detail cover a very wide 
range of sub-disciplines of physics, from atomic physics to equations of state, 
non-ideal chemistry, warm dense matter, particle transport, radiative opac- 
ities, computer simulations of dense matter, computational fluid dynamics, 
and high energy density experiments; all of which are necessary to advance 
the astrophysics of white dwarfs and address the new questions that will 
surely flow from new powerful observatories such as the James Webb Space 
Telescope, the CASTOR space telescope, the Vera Rubin Observatory, the 
SDSS-V, 4MOST, DESI and WEAVE multi-object spectroscopic sky sur- 
veys and the ongoing mining of the data from the Gaia mission. Advances 
in solving the problems we have discussed will have a direct impact in sev- 
eral areas of white dwarf astrophysics. Line profiles and line broadening, 
continuum and molecular opacities, ionization and dissociation equilibria are 
all essential ingredients in modeling atmospheres and for the calculation of 
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model spectra that are directly compared to observations to determine Tyg, 
surface gravity, composition and magnetic field strength. These basic stellar 
parameters inform the tool of cosmochronology to determine the age of stel- 
lar populations and contribute to the timeline of the formation of structures 
in our galaxy. Diffusion, convective mixing, overshoot and thermohaline con- 
vection are all critical to interpret the abundances of the white dwarfs whose 
atmospheres are polluted by infalling planetary solids and to draw conclu- 
sions about the nature of old exoplanetary systems. Crystallization coupled 
with diffusion in the envelope and thermal conductivity affect the cooling 
time scale of old white dwarfs and the determination of white dwarf ages. 
Accurate collision-induced opacities are crucial to understand the sub-class 
of IR-faint white dwarfs about which we know so little. Better models of 
the pressure and magnetic distortion of the Swan bands of C2 and of the 
Cy dissociation in dense He will allow an accurate characterization of the C 
dredge up mechanism that involves both diffusion at the C/He interface and 
convective mixing. Equations of state, in the warm dense matter regime of 
partial ionization of H, He and C affect the extent of convection zones, the 
driving of pulsation modes, and the onset of convective coupling that is a 
critical event in white dwarf cooling. The ionization state of low-Z elements 
is an important factor in inter-diffusion of trace elements, which arcs back to 
the composition of exoplanetary material orbiting white dwarfs. 

It remains a source of fascination that the application of fundamental 
physics principles at the atomic scale, combined with clever models, powerful 
computing resources and experiments on small, laboratory-scale samples can 
lead directly to an understanding of systems as unreachable as stars. In this 
endeavor, we look forward to the contributions of scientists with new ideas 
from related fields of research to further our knowledge of these very small, 
dim, yet numerous stars imbued with cosmic significance. 

We dedicate this review to the memory of our teacher, colleague and friend 
Gilles Fontaine who, animated by a great passion for stellar astrophysics and 
training young scientists, left an indelible legacy in this field. 
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